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1.0 Introduction 


The Personnel Launch System (PLS) Autoland Development study has been performed as 
Task Assignment 8 under NASA contract No. NASf -18762, "Aircraft and Spacecraft 
Guidance and Control Technology". This work was a follow-on effort to that reported in 
reference 7. 

1.1 Background 

A recent NASA effort to provide an Assured Crew Return Capability (ACRC) for the Space 
Station Freedom has resulted in the proposal of several concepts, one of which is the 
Personnel Launch System (PLS) vehicle designed by NASA Langley Research Center. 

The PLS vehicle is to be used to carry up to 10 astronauts to and from the Space Station. 
A requirement for such a vehicle is the capability to perform the entire flight 
automatically, from deorbit burn through reentry to landing. This will provide 
increased reliability and allow for safe return of crews that are either pilotless or in 
which the pilot-qualified crew members are incapacitated. In addition, the vehicle must 
be a fully or partially reusable, and it must provide for emergency crew escape from the 
Space Station and return to Earth and land on a prescribed runway without subjecting 
the passengers to detrimental accelerations. 

The purpose of this task is to develop a practical autoland guidance and control system 
for the PLS vehicle. The autoland system is to be capable of both automatic and manual 
control. This task focuses on the automatic control system and covers two major 
automatic flight phases: the Terminal Area Energy Management (TAEM) phase, and the 
Approach and Landing (A/L) or autoland phase. 

The purpose of the TAEM phase is to transit the vehicle from the high energy state at the 
end of the atmospheric re-entry phase, to the initiation of the approach and landing 
phase with sufficient total energy for a safe landing. The A/L phase initiates at the 
beginning of the turn onto the runway centerline and ends at touchdown. 

The aerodynamic data and vehicle mass properties for the 10 man PLS vehicle were 
provided by NASA Langley Research Center. The aerodynamic data was obtained at Mach 
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0.08 and 0.6, at various angles of attack, sideslip angles and different control surface 
positions. 

The mission requirements and limits are summarized in Table 1-1. These requirements 
are not fixed "hard" requirements, but are intended to serve as a general guideline 
during this preliminary development of the PLS autoland system. For example, the 
requirement that the sink rate be less than 2 feet per second at touchdown is severe and 
may be relaxed, given that the landing gear used on the Shuttle allow sink rates in excess 
of 4 fps. 

The purpose of this study was to develop a control system with the capability of meeting 
requirements, especially tracking the trajectory in the presence of wind turbulence and 
wind shear with both head and side wind conditions. The study also specified that the 
autoland design be verified using a non-linear 6-DOF flight simulation. 

1.2 Technical Approach 

The PLS vehicle was geometrically scaled up in size from CERV (Crew Emergency Rescue 
Vehicle). As a result, the work performed in developing the CERV autoland system (as 
reported in reference 7), was used as a starting point in the design of the PLS autoland 
system. 

The PLS autoland development study focused on developing the guidance and control 
system for performing an automatic landing at a prescribed runway. The study consisted 
of the following five prescribed tasks: 

1 ) Develop TAEM and A/L trajectory profile and algorithms; 

2) Design a multi-input, multi-output (MIMO) system control law; 

3) Develop a non-linear simulation; 

4) Determine autoland guidance and control system feasibility ; 

5) Perform parameter sensitivity studies to evaluate the effect of varying 

mission parameters on vehicle flight performance. 
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Touchdown Requirements 


Requirement 

Definition/Com ment 

a i i5.o Q 

angle-of-attack limit 
( 1 8° hard limit w/ wind) 

T] crab i 5.0£> 

runway heading alignment 
for decrab maneuver 

KEAS*200 knots 

velocity - knots of 
equivalent airspeed 

nominal Zdot= 2.0 fps 
with IQ - lfps 

sink rate 

AY- ±25 ft 

allowable lateral error at 
touchdown 

aimpoinl-2000 ft 

delta-X aimpoint from the 
end of the runway 


Table 1-1 PLS General Mission Requirements/Guideline 
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Our technical approach consisted of developing a MIMO autoland system including pitch, 
yaw, and roll autopilots, lateral steering, altitude and speed controller. The autopilots 
were designed using the integral-LQR control technique. The integral-LQR technique 
consists of using state weighting for transmission zero placement which allows direct 
design in the frequency domain. The resulting design, by using high gain on the integral 
control offers excellent low frequency characteristics with robustness to parameter 
biases. The control weighting is based on decoupling the control surface effectors. The 
design process uses a discrete LQR design technique with computational delay 
compensation. Details of the design synthesis are contained in section 3.0 . 

The resulting MIMO autoland system was verified using a non-linear 6 DOF simulation. 
A block diagram for the simulation is shown in figure 1-1. The CERV 6- DOF simulation 
developed under the previous NASA task (reference 7) was used as a basis for developing 
the PLS 6 DOF simulation. PLS vehicle characteristics, aerodynamic model and the 
digital control system were incorporated into this simulation. 

The stability, control, and performance of the PLS vehicle was verified using this 6 DOF 
simulation. In addition, mission performance was verified for various wind conditions 
including wind shear, wind gusts, and wind turbulence. The 6 DOF simulation was also 
used for parameter sensitivity studies, varying glideslope angles, initial velocity, 
guidance parameters, final flare altitude and speedbrake settings. Details of the 
simulation analysis are described in section 4.0 and 6.0. 


The A/L and TAEM autoland phases were developed separately, and therefore are 
presented separately in sections 3.0 and 5.0 respectively. Section 2.0 describes the 
overall mission, vehicle characteristics and the control techniques.. The A/L system is 
presented first, followed by the TAEM system. 
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Gains scheduled on: 

Mach (M) 
angle-of-attack (a) 
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Angle of attack and sideslip angle 






1.3 Summary of Results 


An autoland system for the PLS vehicle has been developed for the approach and landing 
(A/L) phase, and the Terminal Area Energy Management (TAEM) phase. The feasibility 
of the design of the autoland system has been demonstrated through simulation. 

The baseline control consists of using the Integral-LQR autopilot design with classical 
guidance laws. An overall block diagram is shown in figure 1.-2. Landing reference 
information for the guidance and control is received either from a microwave landing 
system, the global-positioning system (GPS), or both, supplemented with inertial 
reference data. The guidance, navigation and control will use a common set of strap- 
down sensors. The guidance forms angle-of-attack and bank angle commands for the 
autopilot. The modern control autopilot together with fin mixer equations, form the 
control surface actuator commands. The design is an all digital design; sample times and 
computational delays selected are shown in figure 1-2. Details of the autoland system 
are described in section 2.0 . 

A straight-in trajectory profile was developed for the approach and landing phase. The 
trajectory is shown in figure 1-3 and consists of a "double flare” landing profile 
similar to the Shuttle's landing profile. A unique feature of the profile is the pre-flare 
maneuver which is used to decrease the sink rate, modulate the airspeed and reduce the 
flight path angle from a steep glideslope to a shallow glideslope. The double flare landing 
profile is preferred by astronauts/pilots because it presents a less threatening ground 
approach when compared to the single parabolic profile. Starting at the end of the steep 
glideslope, there is approximately 15 seconds to touchdown during which the vehicle 
descends from 300 feet along a flight path angle of -2.5-. The excellent touchdown 
performance obtained using the double flare profile is graphically summarized in figure 
1-4. The sink rate, altitude profile and flight path angle and angle-of-attack during the 
final approach and final flare are plotted in figure 1-4 . This data was obtained from a 
6-DOF simulation of a nominal A/L mission without wind. The sink rate quickly 
dissipates during the final flare, without commanding large angle-of-attack. 
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Autoland System Block Diagram 





Figure 1-3 Summary of PLS Straight-In Trajectory Profile 







Figure 1-4 


Summary 
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The PLS vehicle control configuration and autoland system meets the accuracy 
requirements and desired stability in the presence of wind disturbance. A summary of 
the results in shown in figure 1-5. The touchdown performance with wind turbulence 
was evaluated at five different wind conditions; the results are excellent with all 
requirements being met. 

A TAEM trajectory was developed using a Heading Alignment Cone (HAC) spiral profile 
similar to the Shuttle's TAEM system. The TAEM guidance and control system was 
adapted from the A/L system. This was done to maintain a consistent set of guidance and 
control algorithms and hardware, and to improve transition from the TAEM phase to the 
A/L phase. 

The TAEM trajectory profile developed for this study is shown in figure 1-6. The study 
focused on using a Heading Alignment Cone (HAC). The TAEM spiral HAC trajectory 
algorithms were generated to follow a constantly decreasing turn radius and velocity 
while maintaining a constant bank angle. A 6-DOF simulation was developed to evaluate 
the TAEM performance. The TAEM simulation results with wind turbulence are 
summarized in figure 1-7. These results indicated acceptable TAEM performance. 
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TAEM Wind Input 

Value 

Description 

a x = o Y = a z 

17.3 ft/sec 

wind stnd deviation ; 

L X = Ly= l Z 

1753 ft 

integral scale factor 


TAEM End Condition 

Nominal 

Wind Turbulence 

Description 

AX err 

33 ft 

66 ft 

along track error 

AY err 

-26 ft 

-25 ft 

cross track error 

H 

14620 ft 

14612 ft 

altitude 

11 

= 0 . 8 

= 0 . 2 

heading angle 


Figure 1-7 Summary of TAEM End Conditions With Wind Turbulence 







2.0 PLS System Description 


The following section describes the PLS autoland mission and the control configuration 
and techniques used for the PLS autoland operation. 


2.1 PLS TAEM and Autoland Mission 

The PLS mission is to return astronauts from low Earth orbit to Earth, landing at a 
prescribed runway. The mission consists of de-orbit, re-entry, terminal area energy 
management (TAEM), and finally approach and landing (A/L) on the runway. Only the 
TAEM and the A/L phases of the PLS mission are considered in this study. 

The TAEM trajectory was developed specifically for the PLS vehicle, based on the 
Shuttle's TAEM system. The TAEM flight phases are shown in figure 2-1 . The TAEM 
mission initiates TAEM acquisition after re-entry at approximately 82000 feet and 
Mach 2.5. This phase is performed to acquire a heading which is tangent to the heading 
alignment cone. The heading alignment phase is initiated at approximately 45000 feet 
and Mach 0.95 and performs a constant bank angle maneuver to fly the vehicle around an 
imaginary inverted cone called the Heading Alignment Cone (HAC), until the vehicle is 
aligned with the runway centerline at approximately 15000 ft and Mach 0.6 . At this 
point the A/L trajectory is captured. 

The A/L trajectory was developed for the PLS vehicle, based on the Shuttle's A/L 
trajectory. The A/L flight profile is shown in figure 2-2. The A/L trajectory starts at 
an altitude of approximately 15000 feet and Mach 0.6, with an initial steep glideslope. 
A constant "g" pullup maneuver is performed at an altitude of approximately 2300 feet 
to decrease the sink rate and to capture the final shallow glideslope. The constant "g" 
pullup maneuver is a preflare maneuver that results in a less threatening ground 
approach compared with the standard single parabolic flare maneuver. At approximately 
70 foot altitude, a final flare maneuver is performed to bring the vehicle down to a 
nominal touchdown at 200 knots of equivalent airspeed and a sink rate less than 2 fps. 



• Terminal Area Energy Management (TAEM) Phases: 
1) Acquisition 
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Figure 2-1 TAEM Guidance Phases 
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Figure 2-2 PLS Approach and Landing Trajectory Phases 



2.2 PLS Vehicle Configuration 


The PLS configuration is a relatively low lift-to-drag (max L/D is 3.2 with gear up) 
lifting body configured to perform autoland on one of several commercial landing 
facilities. Figure 2.-3 shows a drawing of the configuration. This configuration is 
similar to the CERV vehicle (see reference 7), but reflects a 15% geometric scaling up 
of the CERV vehicle for increased payload and more passengers. Included in figure 2-3 
is a summary of the vehicle mass properties, comparing PLS mass properties to CERV 
mass properties. 

For control, the vehicle has seven control surfaces as shown in figure 2-4, including 4 
body flaps, 2 wing elevons, and an all-moveable tail for rudder control. The surface 
mixer equations are included in this figure, showing the usage of control surfaces for 
vehicle control. Pitch control can be achieved using either the body flaps or wing 
elevons. Our baseline approach is to use the wing elevons for pitch control and reserve 
the body flaps for roll control and drag modulation. Yaw control is achieved using the 
rudder. Roll control can be achieved using differential body flaps and/or differential 
wing flaps. However, the usage of differential wing flaps for roll control causes large 
yawing moments due to the large wing dihedral. This in turn requires a larger rudder to 
trim the induced yawing moments. This configuration also results in additional drag. 
Therefore, differential body flaps were chosen to control roll. 


2.3 Functional Description of Control System 

The proposed autoland system is an all digital control system consisting of sensors, flight 
processor, interfaces and control surface actuators. The algorithms associated with 
navigation, guidance, and control are contained within the flight processor. A functional 
block diagram for the system is shown in figure 2-5. Navigation uses velocity and 
altitude information from strap-down inertial sensors augmented with altitude 
information from a radar altimeter and inertial measurement updates based on Global 
Positioning System (GPS) data. Vehicle information with respect to the runway is 
received from the Microwave Landing System (MLS). An alternate would be to use GPS 
data instead of MLS data, with a data load for runway altitude and location. 
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MASS PROPERTIES: 


• mwi will 

PLS 1 

CERV 

Crew Size 
Weight (lbs) 
Xcg (%) 

Ixx (slug-ft2) 
lyy (slug-ft2) 
Izz (slug-ft2) 
Lref (ft) 
Aref(ft2) 
Span (ft) 

1 0 

1 9 f 1 00 

55.5% 

7512 

33594 

35644 

28.24 

286.45 

13.89 

8 

1 1,193 

54% 

6345 

13994 

17123 

24.583 

215.61 

12.132 


Figure 2-3 PLS Vehicle Configuration and Mass Properties 
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Figure 2-4 PLS Vehicle Control Surface and Mixing Equations 
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Figure 2-5 Autoland System Block Diagram 






A block diagram for the guidance law is shown in figure 2-6. The longitudinal guidance 
is an altitude/altitude-rate system using feedback information from the Microwave 
Landing System and augmented with navigation information. Output of the longitudinal 
guidance is angle-of-attack command for the pitch autopilot. An angle-of-attack 
command limit is included along with an angle-of-attack trim function based on level 
flight. 


The lateral guidance is a cross track/cross track rate system using information derived 
from the Microwave Landing System. The output of the lateral guidance is bank angle 
command. 


An integral-plus-proportional inertial velocity control loop is also included. Inertial 
velocity feedback information is derived from the navigation system. The speed loop 
commands all four body flaps, splitting them to produce a drag force with a small 
resulting pitching moment depending on the vehicle angle of attack. The pitching moment 
is trimmed by wing elevons. The nominal setting selected was 20 ? during the steep 
descent to account for head winds. 

A block diagram of the autopilots is shown in figure 2-7. The design was developed using 
the integral-LQR technique. The pitch autopilot is an angle-of-attack autopilot with 
integral control and pitch rate damping. Elevator command feedback (SEC) is used to 

compensate for the computer computational delay. The yaw/roll autopilot is a coupled 
design controlling bank angle. Other feedback variables are sideslip angle, yaw rate, 
bank angle, and roll rate. Rudder command feedback (SRC) and aileron command feedback 
(SAC) are also used to compensate for the computer computational delays. A control 
surface mixer is included to change the elevator, speedbrake, and aileron commands to 
individual surface commands. Elevator, speedbrake, rudder, and aileron command 
limiters are also included. A block diagram of the control surface mixer is shown in 
figure 2-8 
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Figure 2-6 Autoland Guidance (Steering Laws) 
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Figure 2-7 Block Diagram of Pitch and Yaw/Roll Autopilots 
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Figure 2-8 Control Surface Mixer Equations 












3.0 Approach and Landing Design Synthesis 


This section discusses the design synthesis for the A/L autopilot, guidance and speed 
controller. The designs are not optimum, but are adequate to conduct studies and to 
assess the PLS autoland system. 


3.1 Autopilot Design 


Linear models were constructed from the aerodynamic derivative data at various trim 
conditions. During approach and landing, the dynamic pressure (Qpres) varies from 
approximately 300 psf at 15,000 feet and Mach 0.6, to nominally 135 psf at sea level 
(200 knots). The design points were chosen at two phases of the flight: first, with 
speedbrake setting at 20 9 during high speed steep glideslope descent, and second, with 
speedbrake setting at 0 9 during lower speed descent. The autopilot gains were designed at 
the following four trim conditions: 

Condition 1 : Q=300 psf, Mach 0.6 a= 8 9 ; Speedbrake=20 9 

Condition 2: Q=270 psf, Mach 0.47 cc= 5 s ; Speedbrake=20 9 

Condition 3: Q=220 psf, Mach 0.41 a= 12 9 ; Speedbrake=0 9 

Condition 4: Q=140 psf, Mach 0.32 a= 16 9 ; Speedbrake=0 9 


The autopilot was synthesized using a modern Integral-Linear Quadratic Regulator (I- 
LQR) design technique. The pitch and yaw/roll autopilots were designed separately, 
since there is little aerodynamic cross coupling between these axes. The following 
sections present a short description of the design methodology, the structure and 
resulting gains of the autopilots. 


3.1.1 Pitch Autopilot 

A derivative and trim program was first used to obtain pitch trim at the four design 
conditions. The derivatives are calculated around a trim condition, using the body or 
wing flaps for trim. The vehicle was trimmed in pitch using the wing flaps; this 
configuration was shown to be desirable, allowing the body flaps to be used for roll 
control and speed control. 



State Space Equation: dX/dt=[A] X + [B] U 
State Vector: X=[ a, q, I a dt]' 

Control Vector: U=[5elev] {wing flap} 



PLANT MATRIX 
[A] 

INPUT 

MATRIX 

[B] 

Eigen- 

values 

Qpress=300 psf 
alpha =8 8 , Mach 0.6 

-.5480 : 

1.0 

: 0 

-.0661 

0 

-8.3636 j 

-.2150 

: 0 

-7.6572 

-.3815 

speedbrake=20 8 

1.0 : 

0 

: 0 

0 

±i 2.8872 

Qpress=270 psf 
alpha =5°, Mach 0.47 
speedbrake=20 8 

-.6235 : 

1.0 

: 0 

-.0752 

0 

-7.4992: 

-.2446 

: 0 

-6.8659 

-.4340 
±i 2.7319 

1.0 : 

0 

: 0 

0 

Qpress=220psf 
alpha=12 8 , Mach 0.41 
speedbrake=0 8 

-.5699 i 

1.0 

: 0 

-.0652 

0 

-5.1877: 

-.3033 

: 0 

-5.048 

-.4366 

H2.2737 

1.0 i 

0 

: 0 

0 

Qpress=140psf 
alpha=16 8 , Mach 0.32 
speedbrake=0 8 

-.4230 \ 

1.0 

j 0 

-.0458 

0 

••••■•• V 

-2.1558! 

-.2456 

: 0 

-2.6754 

-..3343 

±1.4656 

i.o i 

0 

j 0 

0 


Table 3-1 Pitch Linear Models 
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Using the derivative data generated from the trim program, a continuous linear pitch 
model for each condition was developed as shown in table 3-1. The state vector includes 
angle of attack (a), pitch rate (q) and integral of angle of attack (la ) which is included 

to improve on command following and robustness to parameter variation. A block 
diagram of the pitch control loop is shown in figure 3-1 . 

The autopilot was synthesized on the MATRIXx program using the Integral-Linear 
Quadratic Regulator (l-LQR) design technique. This technique helps the designer to form 
a desired state weighting matrix methodically, instead of by trial-and-error as is done 
with the LQR technique. The l-LQR method involves selecting outputs for weighting, 
which produce transmission zeroes of the LQ regulator to control the position to which 
the closed loop poles will migrate. The primary condition imposed by the l-LQR method 
is that the number of output variables must equal the number of input variables to allow 
the specification of the transmission zeroes. For the pitch autopilot, there is one input, 
so one output equation is used for output weighting. 

The equations and development of the l-LQR gains are given in figure 3.-2; the example 
is obtained from the development of the pitch gains for condition 1. First, the 
transmission output matrix [C] is developed from a given set of transmission zeroes. 
Given the cost function J, the state weighting matrix [Rxx] is developed using the 
transmission output matrix [C]. The input weighting matrix [Ruu] is given a nominal 
value of one, and [Rxx] is multiplied by coefficient p, to apply a relative weight between 
[Rxx] and [Ruu]. MATRIXx was used to compute the LQR optimal state-feedback gains 
and the resulting eigenvalues. The coefficient p is varied until desired closed loop 
eigenvalues are obtained. 
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Figure 3-1 Pitch Autopilot Block Diagram 


28 





Using this methodology, the complete continuous gain schedule was developed for each 
condition as summarized in table 3-2. 


Q 

(PSf) 

K« 

Kq 

K », 

Closed Loop Poles 
(continuous) 

300 

-3.890 

-1.188 

-8.174 

-2.717 

-3.70 1±i3. 063 

270 

-3.412 

-1.127 

-7.311 

-2.539 

-3.162±i 3.134 

220 

-4.590 

-1.467 

-8.954 

-2.620 

-2.980± i 2.905 

140 

-6.171 

-2.297 

i 

-9.708 

-2.633 

-2.232± i 2.220 


Table 3-2 Pitch Continuous Autopilot Gains 


The major advantage of using the l-LQR design method over the classical design is that 
the l-LQR design results in larger integral gains, while maintaining larger gain and 
phase margins. This is evident in the integral alpha gains (Kctj) obtained using the I- 
LQR design, as summarized in table 3-2 . 

The discrete pitch autopilot was synthesized using the l-LQR technique. The discrete 
design was aided by using the same LQ weighting matrix developed for the continuous 
case. A 20 msec autopilot sampling rate was used with a 10 msec computational delay as 
shown in the block diagram of figure 3-1. The resulting regulator gains for the discrete 
autopilot are given in table 3-3. 

Discrete gains were also computed for a NASA cockpit simulation study. NASA's 
simulation study for evaluation of the PLS vehicle utilizes the Transportation Systems 
Research Vehicle (TSRV) simulation cockpit and real time simulation facility. This 
simulation facility runs on a 30 msec sampling interval and the computational delay can 
be approximated at 30 msec. As a result of the 30 msec requirement, the discrete pitch 
gains were recalculated at the larger sampling 


29 





• STATE EQUATION: [dX/dt]= [A]* [X] + [B] * [U] 


da/dt 


-0.548 

1.0 

67 


m mm 

a 


”- 0 . 0 6 6~ 

dq/dt 

— 

-8.364 

-0.215 

0. 


q 

+ 

-7.657 

dai/dt 


_ 1.0 

0.0 

0. 


a i 


0 . 


• OUTPUT EQUATIONS TO PLACE TRANSMISSION ZEROES 

Y(s)= fs 2 +2<o$ s + to 2 ) a (s) lei to=3.40 rps; £=0.95 

s a=2co£=6.46; b=co 2 =11.56 

or, Y=(da/dt) + aa+ baj; where: da/dt= A(1,1)*a + A(1,2)*q + A(1,3)*aj 

Y= (A(1 ,1 )+ a)*a + A(1,2)*q + b*aj 
Thus, for condition #1 : 

[Y]= [C] * [X] where: [C]= [A(1,1)+ a A(1,2) b } 

[Y]= [5.912 1. 11.56]* [a q aj] T 


• COST FUNCTION: J= 0 J°° (X T Rxx X + U T Ruu U + U T Rxu X) dt 


. WEIGHTING MATRIX 


j Rxx 

Rxu t 

where; 

Rxx= p [C T C] 

Ruu=[1] 


1 Rxu 

Ruu 


Rxu= [0 0 0] 
Rxx= p[ C^C] = p 

34.95 5.91 
5.91 1.0 

68.34 11.56 

68.34 

11.56 

133.63 


• RESULTING LQR GAINS: 

for p= 0.5 Kr = [-3.89 -1.19 -8.17] 


Eigenvalues= -2.72; -3.70 ±_ i 3.06 


FIGURE 3-2 Example of Integral-LQR Design 
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interval of 30 msec with a computational delay of 30 msec. These gains are also 
summarized in table 3-3. 


Q 

(psf) 

T-samp/ 

T-comp 

Kcc 

Kq 

K «i 

He 

300 

20ms/10ms 

-3.471 

-1.150 

-7.456 

.089 


30ms/30ms 

-3.201 

-1.173 

-7.105 

.266 

270 

20ms/10ms 

-3.088 

-1.099 

-6.757 

.077 

30ms/30ms 

-2.875 

-1.123 

-6.486 

.230 

220 

20ms/10ms 

-4.226 

-1.437 

-8.297 

.074 


30ms/30ms 

-4.019 

-1.473 

-7.978 

! .223 

140 

20ms/10ms 

-5.828 

-2.261 

-9.109 

.062 


30ms/30ms 

-5.680 

-2.318 

-8.816 

.188 


Table 3-3 Pitch Discrete Autopilot Gains 


A second order, 10 Hz actuator model with 0.5 damping ratio was added to the pitch 
autopilot design model. Nichols plots for the 4 conditions with the loop open at the 
actuator are shown in figure 3-3. Plotted are the results from the continuous pitch 
gains, and two discrete conditions: 

-discrete condition 1= 20 msec sampling /1 0 msec computational delay 
-discrete condition 2= 30 msec sampling /30 msec computational delay 

A region of typical stability margin requirements is also shown in the Nichols plots. 
The frequency response results show that there is sufficient gain and phase margin to 
guarantee stability and robustness. However, as expected, the stability margin of the 
discrete case decreases compared to the continuous case, and the stability margin further 
decrease using the larger sampling interval and computational delay. 

In addition, system response to an a step command is plotted for a.q, d/dt(8e) and 8e in 
figures 3-4 for all four flight conditions. The results show a quick time response with 
the actuator elevator rate < 50 deg/sec . 





qp nivo 


qp Nrvo 


o 



qp nivo 


o 



O O O O O O I 
K) <M •— — 04 

I I 


qp nivo 


32 


Figure 3-3 Continuous and Discrete Pitch Nichols Plots 







Figure 3-4 System Response to Alpha Step Command 















3.1.2 Yaw/Roll Autopilot 


A derivative and trim program was used to obtain yaw/roll stability derivatives at four 
flight conditions listed in section 3.1 . The vehicle was trimmed using only the body flaps 
for roll control, and the smaller of two possible all-moveable tail sizes as a rudder for 
yaw control. 

The yaw/roll continuous linear model data for each of the four flight conditions is given 
in table 3-4. The state vector X=[ P r p <}> ] and the input vector U=[ 5 rud 5 ail ]. 

For the PLS design, only the body flaps are used for roll control, and the small all- 
moveable tail as a rudder for yaw control. Previous studies on the CERV vehicle 
(reference 7) showed that using the body flap alone lacked sufficient roll control 
authority to decrab the vehicle with a side wind greater than 15 knots. However, the 
configuration was acceptable for PLS because the requirements for decrab were relaxed 
somewhat by specifying that the heading alignment at touchdown be < 5 s as opposed to 
complete alignment for CERV. Decrab design and analysis was not performed for PLS 
because NASA postponed this for a later study when new data becomes available for 
landing gear and ground effect. 

The yaw/roll autopilot design is based on the CERV design. A block diagram of the 
resulting yaw/roll autopilot is shown in figure 3.-5. The yaw-roll Integral-LQR design 
was done for a state vector that includes sideslip angle, yaw rate, roll rate and roll angle. 
As was done with the pitch autopilot, a continuous autopilot was first designed. The 
resulting gain schedule and eigenvalues are given in table 3-5. 

The same LQ weighting developed in the continuous synthesis was used to obtain a 
corresponding discrete design with computational delay compensation included. Two sets 
of discrete gains were obtained for different sampling rates and computational delay: 

• discrete condition 1: T-samp!e= 20 msec; T-computational delay= 10 msec 

• discrete condition 2: T-sample= 30 msec; T-computational delay= 30 msec 

The corresponding discrete gain schedule is given in table 3-6. 
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***** CONDITION 1 ***** 


V=635 ft/sec; QPRES=300 psf; MACH=0.6; Alpha-5 deg 


-0.1624 

- 1.0000 

0.0873 

0.0000 


0.0267 

-0 .0020 

5.5647 

-0.3082 

0.1254 

0.0000 

B- 

-2.5041 

0.0038 

-71.8552 

1.0638 

-0.9901 

0.0000 


2.4037 

6.0456 

0.0000 

0.0000 

1.0000 

0.0000 


0 .0000 

0 .0000 


Eigenvalues- 0.; .1648; -.8127-3 . 487i ? -. 8127+3 . 487i 
***** CONDITION 2 ***** 


V-500 ft/sec; QPRES-270 psf; MACH-0.47; Alpha-5 deg 



-0.1856 -1.0000 0.0873 0.0000 


0.0306 -0.0023 

A= 

5.0082 -0.3523 0.1433 0.0000 

B« 

-2.2537 0.0035 


-64.6697 1.2159 -1.1316 0.0000 


2.1633 5.4410 


0.0000 0.0000 1.0000 0.0000 


0.0000 0.0000 


Eigenvalues- 0.; .1831; -0.9263 - 3.2607i; -0.9263 + 3.2607i 


***** CONDITION 3 ***** 


V-435 ft/sec; QPRES-220 psf; MACH-0.41; Alpha-10 deg 



-0.1738 -1.0000 0.1745 0.0000 


0.0335 -0.0042 

A- 

5.6286 -0.3612 0.0859 0.0000 

B- 

-1.9911 0.0619 


-52.6938 1.5630 -0.7406 0.0000 


2 . 5305 4 . 1664 


0.0000 0.0000 1.0000 0.0000 


0.0000 0.0000 


Eigenvalues- 0.; -.0967; -0.5895 + 3.8325i; -0.5895 - 3.8325i 


***** CONDITION 4 ***** 


V-340 ft/sec; QPRES-140 psf; MACH-0.32; Alpha-16 deg 


-0.1415 

- 1.0000 

0.2792 

0.0000 


0.0326 

-0.0048 

2.7759 

-0.2944 

0.0533 

0.0000 

B- 

-1.3665 

0.0609 

-33.5324 

0.8666 

-0.7560 

0.0000 


2.0437 

2.7873 

0.0000 

0.0000 

1.0000 

0.0000 


0.0000 

0.0000 

— 









Eigenvalues —3.9417 - 1.3138i; -3.9417 + 1.3138i 
-3.8689 + 3 . 7090 i ; -3.8689 - 3.7090i 

. lT 

JC= C*\X + Ch]UL pOfykrj, 4>u^9l 

VL- dVu4 6*g), £axl(/*$)\ T 

TABLE 3-4 YAW/ROLL LINEAR MODEL DATA 
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Figure 3-5 PLS Yaw/Roll Autopilot Block Diagram 
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K- RUDDER CONTINUOUS GAINS 



Condition 

# K-beta 

K-r 

K-p 

K-phi 

i 

11 . 5602 

-3.1173 

-0.2505 

-0.8380 

2 

11.6715 

-3.2571 

-0.2791 

-0.9199 

3 

10.9234 

-3.4077 

-0.2393 

-1.2800 

4 

12 . 5104 

-4.3575 

-0.0294 

-1.9237 

K- AILERON CONTINUOUS GAINS 



Condition 

# K-beta 

K-r 

K-p 

K-phi 

1 

-14.2882 

1.7535 

1.8251 

5.9263 

2 

-14.1947 

1.8542 

1.9171 

6.2254 

3 

-15.5621 

2.5746 

2.2310 

6.6544 

4 

-16.3209 

3.3596 

2.8141 

8.8851 


TABLE 3-5 YAW/ROLL CONTINUOUS GAIN SCHEDULE 


K-RUDDER DISCRETE 

GAINS (Tsamp- 

20 msec; 

Tcomp-10 

msec ) 


Condition # K-beta 

K-r 

K-p 

K-phi 

K-dar 

K-drr 

1 9.9224 

-3.0850 

-0.1863 

-0.6318 

0.1462 

-0 . 0232 

2 10.1188 

-3.2275 

-0.2137 

-0.7026 

0.1383 

-0 .0239 

3 9.3022 

-3.3565 

-0.1649 

-1.0523 

0.1281 

-0 .0189 

4 10.5550 

-4.2032 

0.0104 

-1.6707 

0.1192 

-0.0061 

K -AILERON DISCRETE 

GAINS (Tsamp 

-20 msec; 

Tcomp-10 

msec ) 


Condition # K-beta 

K-r 

K-p 

K-phi 

K-daa 

K-dra 

1 -15.7416 

2.0054 

1.6420 

5.0488 

-0.0224 

0.1964 

2 -15.5594 

2.1061 

1.7343 

5.3467 

-0.0225 

0 .1869 

3 -16.6148 

2.8430 

2.0109 

5.7949 

-0.0156 

0.1707 

4 -16.8250 

3.5638 

2.4574 

7.5560 

-0.0032 

0 .1419 


K-RUDDER 

DISCRETE 

GAINS ( Tsamp= 

30 msec; 

Tcomp=30 

msec ) 


Condition 

# K-beta 

K-r 

K-p 

K-phi 

K-dar 

K-drr 

1 

9.8524 

-3.1732 

-0.1882 

-0.6400 

0.2224 

-0.0351 

2 

10.0612 

-3.3167 

-0.2162 

-0.7104 

0.2102 

-0 .0362 

3 

9.1642 

-3.4356 

-0.1633 

-1.0595 

0.1943 

-0 .0285 

4 

10.4016 

-4.2920 

0.0184 

-1.6763 

0 .1808 

-0 .0089 

K-AILERON 

DISCRETE 

GAINS (Tsamp 

=•30 msec; 

Tcomp=30 

msec ) 


Condition 

$ K-beta 

K-r 

K-p 

K-phi 

K-daa 

K-dra 

1 

-16.7603 

2.1770 

1.6568 

5.0204 

-0 .0398 

0 . 2954 

2 

-16.5259 

2.2780 

1.7498 

5.3198 

-0 .0394 

0.2812 

3 

-17.4681 

3.0330 

2.0202 

5.7737 

-0 . 0292 

0.2566 

4 

-17.5199 

3.7460 

2 .4628 

7.5383 

-0 . 0087 

0.2132 


TABLE 3-6 YAW/ROLL DISCRETE GAIN SCHEDULE 
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To evaluate the stability margins, Nichols plots with the control loops opened 3t each 
3Ctu3tor were generated for the 4 flight conditions, and ore plotted in figures 3-6 
through 3.-9. These plots include the actuator dynamics ss described for the pitch 
autopilot analysis. Each figure contains the open yaw loop Nichols plot with the roll loop 
closed, snd the open roll loop Nichols plot with the yaw loop closed. Plotted are the 
results from the continuous yaw/roll gain, snd for the two discrete conditions. As 
expected, the stability margins of the discrete esse decrease compered to the continuous 
esse, snd the margins decrease further using the larger ssmpling time 3nd Isrger 
computstionsl delsy. However, all stability margins are sufficient and do not violate the 
design margins shown on the plots. 

In addition to the frequency analysis of the yaw/roll autopilot, a time response study was 
performed. The system response to a roll step input was evaluated at each of the four 
flight conditions. The results are plotted in figure 3.-10. Only 1.5 seconds is needed for 
the roll response to equal the roll command. 
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Figure 3-9 Yaw/Roll Nichols Plot for Flight Condition #4 
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3.2 A/L Longitudinal Guidance and Control System 


The A/L longitudinal guidance and control system is made up of the trajectory (or flight 
profile) generator, the altitude control module and the speed control module, as shown in 
the signal flow diagram of figure 2-6. The trajectory generator determines the desired 
altitude and descent rate from the vehicle current position relative to the runway. 
Trajectory algorithms were developed to calculate the PLS flight profile. The guidance 
system operates in a 60 msec time frame. 

The along track error is indirectly controlled by the altitude control module and the 
speed controller. The altitude controller translates altitude and sink rate commands 
from the trajectory generator into angle-of-attack commands. The speed controller is 
covered in section 3.4. 


3.2.1 A/L Trajectory Generator 

The A/L trajectory generator generates straight-in trajectory commands. The PLS 
straight-in trajectory uses a Shuttle-like landing profile. The trajectory and algorithm 
variables are defined in figure 3-11. The trajectory profile can be divided into 5 
sections: 

1- steep glideslope (HO to HI, glide slope=-26 9 ); 

2- pull-up (or pre-flare) maneuver (constant g circle;H1 to H2C); 

3- exponential transition (shallow glideslope capture , H2C to H2) 

4- shallow glideslope (H2 to H3, glide slope=-2.5 Q ) 

5- final flare and decrab (H3 to touchdown). 

The trajectory algorithms calculating the trajectory variables is given in appendix G. 
Figure 3-11 also defines preliminary trajectory parameters. These parameters define 
a nominal PLS mission. They were determined after repeated 6-DOF simulations were 
performed, varying the input parameters until satisfactory landing performance was 
obtained (see section 4.1.1). The resulting mission initial conditions that define the 
nominal PLS A/L straight-in trajectory are summarized in Table 3-6. 
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Figure 3-11 PLS Straight-In Trajectory Profile 



Input Parameter 

Value 

Description 

GAMO 

-26® 

steep qlideslope angle 

GAM1 

-2.5® 

shallow qlideslope angle 

WFLAR1 

0.220 

1st flare rate (rps) 

WFLAR2 

0.280 

2nd flare rate (rps) 

ZFLARE 

1.0 

flare damping ratio zeta 

DELS 

7000 


DELH 

2000 

1st flare (constant g pullup) delta range 

HI 

2300 

start of 1st flare altitude 

H2C 

300 

end of 1st flare altitude 

H2 

250 

shallow qlideslope capture altitude 

H3 

70 

start of 2nd (final) flare altitude 

HTDC 

- 1 0 

commanded touchdown altitude 

ACIR 

0.5 

q pullup durinq 1st flare 


Table 3.-6 Mission Input Parameters For A/L Nominal Trajectory 

Each of the 5 trajectory profile sections will be explained in greater detail. The first 
A/L phase starts with a steep glideslope which is strongly dependent on the lift-to-drag 
(UD) characteristics of the vehicle. Due to the low L/D ratio of the PLS vehicle, a glide 
path angle of 24® or higher is necessary to maintain enough total energy to carry the 
vehicle through to touchdown. For the final nominal trajectory, a 26® glide path angle 
was chosen. This allows for a nominal speedbrake deflection of 20® during the steep 
glideslope descent phase necessary to control speed in the presence of wind. 

The second A/L phase is the preflare (or pullup) maneuver consisting of a constant "g" 
pullup maneuver. This maneuver is performed to decrease the sink rate and to capture 
the final shallow glideslope. This preflare maneuver is prefered by astronauts/pilots 
because it results in a less threatening ground approach compared with the standard 
single parabolic flare maneuver. The trajectory generator commands a constant "g” 
circular profile defined by a constant radius. The radius length results in a M g” load 
factor, and is determined from the desired circular delta altitude (DELH) and delta range 
(DELS) as shown in figure 3-11. A 0.5 M g” load was chosen for the PLS pullup. 




































The third phase of the trajectory is the exponential transition phase. A second order 
exponential decay filter is used to generate an exponential transition profile to capture 
the shallow glideslope angle. The algorithm calculates an altitude acceleration command. 
The acceleration command decays to zero as the commanded altitude rate and position 
reaches the desired final value. 

Kinetic energy is quickly lost as the vehicle remains in the preflare (constant g pullup) 
maneuver. To minimize the loss of energy during this maneuver, simulation results 
showed that the final glideslope angle at the end of the preflare should be as large as 
feasible. However, it is also desirable for the final glideslope angle at the end of the 
preflare be equal to the shallow glideslope angle; this minimizes the transients that 
occur during the exponential transition . As a result, it was established that the final 
flight path angle at the end of the preflare should be slightly greater than the shallow 
glideslope angle. By varying the DELH and DELS parameters, the flight path angle at the 
end of the preflare can be controlled. As a result, DELH and DELS were adjusted to 
obtain a glideslope angle of approximately 5 s , which is slightly larger than the shallow 
glideslope of 2.5 9 . 

The shallow glideslope phase is the final low altitude approach to the runway. To 
maintain sufficient kinetic energy, a flight path angle of -2.5 Q is required. 

The final flare profile is similar to the exponential transition flare, and is used to 
reduce the sink rate at touchdown to an acceptable level. The most sensitive parameters 
affecting the touchdown conditions are the final flare rate (WFLAR2), the altitude at 
which to start the final flare (H3) and the commanded touchdown altitude (HTDC). The 
final flare should start at approximately 70 feet to allow the autoland system to decrease 
the sink rate prior to touchdown. A lower H3 altitude results in a higher landing speed 
but at the expense of a larger sink rate. On the other hand, a higher H3 results in 
lowering the sink rate but at the expense of an unacceptably low landing speed. 
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3.2.2 A/L Altitude Control Module 


The function of the altitude control module is to translate the altitude and sink rate 
commands from the trajectory generator into appropriate angle-of-attack signals for 
the flight control. A block diagram of the longitudinal guidance loop is shown in figure 
3-12. The guidance loop is designed classically, and consists of an altitude and altitude 
rate feedback. The gain KH associated with the altitude error forms a zero in the root 
locus plot. Gain KAC is the guidance loop gain. An altitude step response analysis was 
performed using the 6-DOF simulation to evaluate the altitude gains. The nominal value 
for KAC is 0.3; the KH gain is a function of altitude as shown in figure 3-12. 

Trimmed angle-of-attack generated from the lift curve as a function of dynamic 
pressure is summed with the altitude steering command to form a total angle-of-attack 
command. Feed forward loops are used instead of an integrator in the altitude loop, and 
as a result, an altitude error bias is observed when the calculated trimmed angle-of- 
attack is not exactly matched with the real trim angle. However, this error bias is small 
and does not affect the landing performance. An altitude error integrator can reduce this 
bias, but it may pose stability problems in the steering loop due to the presence of an 
additional pole at the origin. 



Figure 3-12 Longitudinal Guidance Design 
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3.3 A/L Lateral Guidance and Control System 


The approach and landing lateral guidance system tracks the runway centerline during 
the descent. The system uses a roll-to-steer method by generating bank angle commands 
to null the cross track error and cross track rate with respect to the runway centerline. 
A block diagram of the lateral guidance loop is shown in figure 3-13. The feedback 
signals are the lateral displacement (Y) from the centerline and the lateral velocity 
normal to the fixed center line. The KY and KPHI gains were varied to evaluate the 

response to a lateral step input. As a result of this study, the best results were obtained 
with KY=0.4 and KPHI=2.0 . 





Figure 3-13 A/L Lateral Guidance Block Diagram 


3.4 A/L Speed Control 


PLS uses a speed controller method similar to the Space Shuttle's. The vehicle is 
controlled to follow a predetermined speed profile. The purpose of the speed controller 
during the A/L phase is to deliver the vehicle to a proper inertial speed at the start of 
the first flare. The desired nominal speed at the initiation of the first flare is 
approximately 350 knots of equivalent airspeed (KEAS). 

The speed controller is most effective during the steep glideslope descent phase. During 
this phase, the speedbrake nominal setting is 20?. This setting allows for speedbrake 
deflections of ±20 e to account for head and tail winds. The speedbrake retracts to a zero 
bias setting at the initiation of the first flare which occurs at an altitude of 2300 feet. 
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This results in a minimum drag configuration which is needed as the total energy quickly 
dissipates during the remainder of the flight. 

A block diagram of the speed controller is shown in figure 3.-14 . The integral gain 
(K51) is used to reduce the regulator error (velocity error). The integrator stops 
integrating when the speedbrake position is limited. The upper limit is set by defining 
the input parameter SPBRUL (=30 s nominal). The lower limit defaults to zero. The 
speedbrake bias (5sb-bias) * s set steep glideslope descent phase at 20-, and for 

the remainder of the flight the bias is set to 0 9 . 

The 6-DOF simulation was used to evaluate the regulator gains and the integral gain 
values. The best results were obtained with the following gains: 

• Ks = -1.2 

• Ksi = 0.4 

Wind conditions will affect the speedbrake performance. The wind uncertainty is 
corrected by the speed controller primarily during the steep glideslope descent . With a 
head wind, the speed controller retracts the speedbrake to reduce drag and maintain a 
constant KEAS. The speed controller performance in the presence of a maximum 
expected head wind of 22 knots was evaluated. A plot of the wind, speedbrake command 
and the velocity error is shown in figure 3-15 . These data were obtained during the 
steep glideslope descent phase. The results show the velocity error is well regulated, 
despite the speedbrake being limited at zero deflection during the last few seconds. The 
speed brake performance is covered in greater detail in section 4.5 . 



Stop integrating when limited 



v 

Figure 3-14 Block Diagram of A/L Speed Controller 
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Figure 3-15 Speed Controller Performance (w/22 knot head wind) 
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4.0 Approach and Landing 6 DOF Simulation Analysis 


The 6-DOF simulation was used to evaluate the performance of the A/L design. First, a 
parameter sensitivity study was performed to evaluate the effect of mission input 
parameters on vehicle touchdown performance. Then, different wind conditions were 
added to the simulation to evaluate the effect of wind on vehicle performance. Finally, 
the speed controller performance was analyzed. 


4.1 Parameter Sensitivity Study 

The 6-DOF simulation was used to evaluate the effect of mission input parameters on 
final touchdown performance. The purpose of this study is to determine what factors 
have the most influence on landing performance, that is, what is the sensitivity of a 
given parameter to landing performance. The results aid in determining nominal A/L 
mission input conditions and also give general trends that can be used in designing a 
mission. For example, increasing the shallow glideslope angle results in an increase of 
the touchdown velocity and sink rate, and decreases angle-of-attack. Therefore, if a 
higher vehicle velocity is desired to provide a larger energy margin in the presence of 
wind uncertainty near touchdown, and larger sink rates are acceptable, the shallow 
glideslope angle may be increased. 

A test matrix was developed to perform repeated 6-DOF simulations without wind, 
varying the following A/L mission input parameters: 


• GAM0= 

steep glideslope angle 

• GAM1 = 

shallow glideslope angle 

• VT0= 

initial true airspeed 

• V1 = 

velocity at start of pullup maneuver 

II 

CM 

X 

• 

shallow glideslope capture altitude 

• H3= 

start of 2nd (final) flare altitude 

• HTDC= 

commanded touchdown altitude 

• WFLAR2= 

2nd flare command filter frequency (rps) 



The following states are used to evaluate touchdown performance: 

• Zdot= Sink rate (fps) 

• a = angle-of-attack (deg) 

KEAS= velocity in knots of equivalent airspeed 
Qatm= Q-atmosphere, or dynamic pressure (psf) 


The results are plotted in figures 4-1 through 4-9. These simulations were run 

without wind. Each plot contains the resulting touchdown conditions for a given mission 

input parameter. The mission general guidelines/requirements at touchdown given in 

able 1 1 serve as a tool in evaluating the landing performance. These mission 

guidelines, applied to landing in the absence of wind, results in the following design 
criteria : 

• Zdot < 2 ft/sec 

• a < 1 5 degrees 

• KEAS > 200 knots 


The effect of the steep glideslope angle (GAMO) on touchdown performance is shown in 
figure 4-1. In general, the increase in the glideslope angle results in better landing 
performance. Note that GAMO-22’ yields an unacceptably high sink rate of 16.6 fps 
The results of this test show that GAM0<-24’, the optimum being GAM0=-26’. 

The effect of the shallow glideslope angle (GAM1) on touchdown performance is shown in 
figure 4-2. The touchdown performance is very sensitive to this parameter. Too small 
an angle (GAM1=-1.5’) results in the vehicle running out of speed and dynamic 
pressure causing a large angle-of-attack and an unacceptably high sink rate at 
touchdown. Essentially, the aircraft loses lift and -falls out of the sky". At the other end 
of the spectrum. (GAM1 « -3.0’), a higher glideslope angle results in sufficient speed 
and dynamic pressure, but also yields a large sink rate (Zdot=3.60 fps). The best 
landing performance meeting all of the design criteria is obtained with GAMl=-2.5’ . 

The touchdown performance is not very sensitive to the change in altitude at initiation of 
the shallow glideslope (H2) as can be seen in figure 4-3. It is desirable to initiate this 
at an altitude high enough to reduce the sinkrate to an acceptable value prior to the final 
flare. On the other hand, too high an altitude may cause unacceptable energy loss and 
stall before touchdown. As a result, a nominal value of H2=250 feet is chosen. 







The sink rate at touchdown is very sensitive to the change in altitude at initiation of the 
fmal flare (H3) as shown in figure 4-4. In order to maintain a low sink rate at 
touchdown, H3 was chosen to be greater or equal to 70 feet; this allows for a sufficient 
altitude for the autoland system to decrease the sink rate. However, increasing H3 has 

the negative effect of lowering the landing speed. Therefore, the nominal value of H3 is 
set to 70 feet. 


Figures 4-5 and 4-6 shows the effect of additional touchdown parameters on landing 
performance. The commanded touchdown altitude (HTDC) has a dramatic affect on the 
touchdown performance, especially on the sink rate. Only one of the four specific 
conditions investigated satisfies all the design criteria, at HTDC=-10 feet. HTDC=-5 
feet yields better sink rate performance but at the expense of a lower KEAS. The final 

flare rate (WFLAR2) also affects the touchdown performance as shown in figure 4-6 
The best choice is for WFLAR2=0.28 rad/sec. 

The next two plots show an interesting contrast between using the speed controller 
versus not using ihe controller, as a function of initial velocity (VTO). As could be 
expected, wilh the speed controller turned oil, the touchdown end condilions are a linear 
lunchon ot the initial velocity as shown in figure 4-7. A 50 ll/sec increase in initial 
velocity yields a corresponding 6 fl/sec increase in final touchdown velocity a linear 
increase in sink rate, and a linear decrease in angle-of-attat*. However, with the speed 
controller turned on. the touchdown condilions remain fairly constant despite a change in 

. initial velocity as shown in figure 4-8 . This result shows the eltectiveness of the speed 
controller in taking out velocity deviations. 

The results of the speed controller comparison test indicate that in order to affect the 
touchdown conditions, the speed controller setting must be changed. The speed controller 
■s most effective during the steep glideslope descent phase, during which a constant KEAS 
velocity (VI) is commanded up to the initiation of the pullup maneuver. A test was 
performed varying VI from 310 to 370 KEAS and the results are plotted in figure 4-9 . 
The results show that for an acceptable touchdown, Vi >330 KEAS should be used. It is 
desirable to maintain a VI low enough to avoid reaching the maximum allowable dynamic 
pressure when wind turbulence or gusts occur. This criteria yields the selection of 
VI 330 KEAS. However, wind disturbance may decrease the touchdown speed and 
require a higher VI. As a result of wind conditions tests. VI is chosen to be V1=350 
KEAS. This case is presented in section 4.3.3 and 4.3.4. 
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4.2 Landing Performance Without Wind 


As a result of the parameter sensitivity study, the mission input parameters that result 
in an acceptable touchdown were determined for the case without wind. These 
parameters, summarized below in table 4-1, form a database that defines the nominal 
PLS autoland mission parameters . These parameters were tested using the 6-DOF 
simulation to evaluate the touchdown performance with the no wind condition. The 
results of the final touchdown conditions are summarized in table 4-1. 

The complete simulation results are given in appendix D. However, a few plots are given 
in figures 4-10 through 4-12 to show general autoland landing performance. Figure 
4-10 shows the velocity and dynamic pressure dynamics as a function of range (s). The 
constant KEAS velocity of 330 knots is clearly seen during the steep glideslope descent 
phase. The altitude profile versus range is shown in figure 4-11. 

Figure 4-12 contains the altitude, sink rate, angle-of-attack, flight path angle and 
normal acceleration dynamics as a function of time, during the final approach and flare 
phases. The results show an excellent sink rate (Zdot) performance. The final flare 
initiates with Zdot=16 fps at approximately 60 foot altitude, and decays such that at an 
altitude of 20 feet, the sink rate is 10 ft/sec. In addition, the angle-of-attack profile 
during the entire approach and landing is benign, requiring approximately 14® at 
touchdown. 

The simulation without wind established the nominal autoland performance of the PLS 
vehicle. The addition of wind turbulence and discrete wind gust is the final test of the 
autoland performance. This is covered in the next sections. 
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cross-track error (feet) 


Table 4-1 Touchdown Performance For A/L Nominal Trajectory 
Without Wind 






















} . AAUVV. A.VVVV* . » 

S tfO 

Figure 4-10 Velocity and Dynami 
A/L Simulation - I 









■ ^11! 

--i -4 1 I 

i » • ' 

-4 1 1 I 

I I I I 

- + -+.- + -1 

C&Jszc)^j£ l i 

-+-+-+-H 

° a + 3 l|- + -M 

:~t t T 1 

g 4^44-i 

hir 1 ' _ + " ^ + + .' Nr^to) 


T+> 


" I 

h«" 


+- 


• + - + --H- + - + 

V_v Ofes) 'I I ' ' 

1 4 - + -.4- 

63. 

sho-llow ^ j 

^Ude slope ' 


Fiaure 4-12 Touchdown Time History for Nominal 
A/L Simulation - NO Winds 


64 



4.3 Landing performance in wind turbulence 


The wind turbulence model is extracted from the work of Holley and Bryson (reference 
4), in the article "Wind Modelling and Lateral Aircraft Control for Automatic Landing", 
January 1975, Standford University Department of Aeronautics and Astronautics Report 
(SUDDAR) No. 489. The model approximates the major characteristics of the Dryden 
and Von Karman wind models. A description of the wind model is included in Appendix C. 

The 6-DOF simulation contains the GRAM atmospheric database with 16 different 
atmospheres. This model was used to generate a steady wind profile. One of the 16 
atmospheres was chosen as the nominal wind profile to be used for evaluating the PLS 
performance with wind. This profile was chosen because it simulated a strong 
longitudinal wind useful for evaluating longitudinal performance. The PLS nominal A/L 
flight trajectory was flown through this steady wind profile and was augmented with the 
wind turbulence data. The final wind profile is shown in figure 4.-13. The wind 
turbulence has a Gaussian distribution of la equal to 5 knots for the x and y axes, and 3 
knots for the z axis. The integral scale length is Lu=69 ft, Lv=36 ft, and Lw=16 ft. 

The effect of the wind turbulence on landing performance is summarized in table 4-3. 
The results show acceptable landing performance with all requirements being met. The 
sink rate of 1.80 ft/sec is excellent. The AX and AY values represent the deviation from 
the nominal runway touchdown without wind, and falls within the required landing 
"footprint". Time history plots are given in Appendix E . 


Touchdown Condition 

Value 

Description 

Zdot 

1 .80 

sink rate (fps) 

KEAS 

180. 

knots of equivalent airspeed 

Vtot 

304. 

ground velocity (fps) 

Qatmos 

1 17 

dynamic pressure (psf) 

a 

17.0 

angle-of-attack (deg) 

ncrab 

= 0. 

crab angle(deg) 

AY 

= o. 

cross-track error (feet) 

AX 

-11 6. 

along-track error (feet) 


Table 4-3 Touchdown Performance With Wind Turbulence 












4.4 Landing Performance in Wind Shear 


A wind shear of magnitude 8 knots/100 ft altitude was added to the wind to obtain a wind 

shear profile at landing. The wind shear was added first as head wind, and then as a cross 
wind. 

Three wind conditions were evaluated. Two head wind profiles were developed as shown 
in figure 4-14. The first profile contains a constant 30 knot head wind profile at 
landing, and second profile contains a 30 knot head wind with a decreasing final shear. 
With both cases, the head wind increases from 0 at an altitude of 4925 feet to 30 knots 
at an altitude 2300 feet, the altitude at which the pullup maneuver begins. The 
increasing headwind at this part of the flight presents a worst case scenario in terms of 
speed control because once the pullup maneuver begins, the speed profile is dominated by 
high deceleration and changing wind conditions. However, this wind would increase the 
airspeed which is not necessarily worst case in terms of maintaining sufficient airspeed 
for the landing. The second wind condition is a worst case condition in terms of landing 
because of the decreasing headwind starting at 375 feet. This will cause a decreasing 
airspeed which could result in unacceptably low touchdown airspeed. The third wind 
shear condition is a cross wind profile. This profile contains an increasing cross wind 

shear of 8 knots/100 ft altitude to a maximum of 22 knots at landing, as shown in figure 
4-14. 

A summary of the touchdown conditions is given in table 4-3 for all three wind shear 
conditions. 

As expected, excellent landing performance was obtained with the constant 30 knot head 
wind (condition 1). The addition of the final shear greatly affected the angle-of-attack. 
This wind shear is similar to having a sudden tail wind. The sudden decrease in total 
velocity due to the wind shear, causes a decrease in lift. As a result, a high angle-of- 
attack is commanded. Despite this, the touchdown requirements were still met. 




Touchdown Results 


Parameter 

Cond. 1 

Cond. 2 

Cond. 3 

Description 

Zdot 

1.84 

1.90 

2.0 


KEAS 

205. 

180 

211 

knots of equiv.airspeed 

Vtot 

346. 

110 

359 

qround velocity (fps) 

Qatmos 

142 

304 

152 

dynamic pressure (psf) 

a 

13.9 

18.0 

13.0 

angle-of-attack (deg) 

qcrab 

= 0. 

= o. 

6.0 

crab angle(deg) 

AY 

=: 0. 

II 

o 

8.0 

cross-track error (ft) 

AX 

215. 

206 

5.0 

along-track error (ft) 


Table 4-3 Touchdown Performance With Wind Shear 


A decrab maneuver was not performed at touchdown. Despite this, the cross wind shear 
of 22 knots resulted in a crab angle of only 6 degrees. 































ctf 




69 


Figure 4-14 Wind Shear Profile 



4.5 Speed Control Performance 


To evaluate the speed control performance in the presence of head winds, simulation 
runs were executed with the speed controller OFF and ON . The results are shown in 
figure 4.-15. The resulting velocity error is plotted for two conditions of speed 
controller off and on. The vehicle inertial velocity is controlled to within 2 fps of the 
velocity command. The wind profile is superimposed on the plot. Note the additional 
discrete gust profile which occurs at the end of the mission when the speed controller is 
disabled. The resulting speedbrake command is shown in figure 4-16 with the wind 
profile included. 


The speed control performance is further proven in the summary of the touchdown 
condition given in table 4-4. The nominal simulation without wind is used as a measure 
of speed controller performance. With head winds shown in figure 4-15, the speed 
controller was turned OFF and then ON. With the speed controller OFF, the touchdown 
speed (KEAS) is unacceptably lowered by 1 1 knots, and the touchdown footprint (AX) is 

dispersed 54 feet from the nominal. However, with the speed controller ON, the 
touchdown sink rate (ZD) and the velocity (KEAS) are identical with the nominal no 
wind case; the touchdown footprint is only dispersed by 4 feet. 


Simulation 

Case 

Speed 

Control 

Wind 

Zdot 

fps 

KEAS 

knots 

Alpha 

deg 

AX-touchdown 

feet 

Nominal Case 

CN 

OFF 

1 .87 

196 

15.4 

N.A. 

ControllerOFF 

OFT 

CN 

1.85 

185 

17.1 

-54. 

Controller ON 

CN 

CN 

1.87 

196 

15.2 

-4. 


Table 4-4 Summary of Touchdown Conditions With Speed 
Controller OFF Versus ON 
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Figure 4-15 Velocity Error With Speed Controller ON Versus OFF 
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Figure 4-16 



5.0 TAEM Design Synthesis 


Our approach for the Terminal Area Energy Management (TAEM) design synthesis is to 
integrate the TAEM mission together with the guidance and control laws used on the 
approach and landing phase. The guidance gains have been developed consistent with the 
TAEM flight conditions and maneuver requirements. Minor changes are incorporated in 
the autopilot gains to assure stability at the TAEM flight conditions. An overview block 
diagram for the integrated TAEM mission generator and guidance law is shown in figure 
5-1. The TAEM trajectory generator forms commands for the longitudinal guidance, 
speed controller.and the lateral guidance. 


5.1 TAEM Guidance and Control 

The TAEM longitudinal guidance system is made up of the trajectory generator and 
altitude controller. The trajectory generator generates a spiral trajectory command. 
The altitude controller translates altitude and sink rate commands from the trajectory 
generator into angle-of-attack commands. 

The following section describes the trajectory generator and longitudinal guidance and 
control. 
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Figure 5-1 Simplied Block Diagram of TAEM Guidance and Control System 






5.1.1 TAEM Trajectory Generator 


Our approach in the PLS trajectory generator is to adapt the Shuttle's TAEM system. The 
Shuttle's TAEM consists of the following 4 guidance phases : 

• Phase 0: S-Tum Phase is an optional phase for excess energy condition. 

Vehicle turns away from the nominal approach to increase ground 
track range. 

• Phase 1 : Acquisition Phase is performed to acquire a heading which is tangent 

to the heading alignment cone (HAC). 

• Phase 2: Heading Alignment Phase rolls the vehicle to follow a spiral 

trajectory along the heading alignment cone. 

• Phase 3: Prefinal Phase converges trajectory to the final approach plane until 

A/L phase is captured within a set criteria. 


Guidance phases one through three are graphically illustrated in figure 5-2 . TAEM 
acquisition phase initiates at Mach 2.5 and 82000 ft altitude. The heading alignment 
phase initiates at Machd.O at a nominal altitude of 35,000 feet and terminates at the 
acquisition of the A/L phase at Mach 0.6 and 15,000 ft altitude. 


In this report, the design analysis is included for phases 2 and 3 (heading alignment and 
prefinal phases). The two initial TAEM phases are left to future analysis. 
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• Terminal Area Energy Management (TAEM) Phases: 
1) Acquisition 
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Figure 5-2 TAEM Guidance Phases 



Figure 5-3 shows the HAC configuration in greater detail. The heading alignment phase 
maneuvers the vehicle onto final approach by flying around the HAC. The HAC is an 
imaginary inverted cone tangent to the runway centerline at approximately 15000 ft 
altitude and 45000 ft from the runway. Way Point 1 (WP1) is the target entry point to 
the HAC; this is the desired vehicle location to initiate the HAC turn. The A/L nominal 
entry point (NEP) is the target entry point to the runway centerline; this initiates the 

A/L capture. 

The prefinal phase simply provides a smooth transition from TAEM to A/L . In prefinal, 
the guidance turns the vehicle heading towards WP2 (Way Point 2) and activates the A/L 
guidance once all capture criteria are met. 

For the Shuttle vehicle, there are 4 optional TAEM targeting (OTT) profiles; each OTT 
results in a different HAC position based on the entry point and on vehicle location. 
Briefly, the option is based on either a straight-in (nearest HAC) or overhead (far HAC) 
approach, and on selecting a nominal or a minimum entry point, as shown in Figure 5-4. 
The minimum entry point (MEP) HAC is selected when a low energy condition exists. In 
addition, the spiral radius during the HAC maneuver can also be shrunk to reduce the 
range-to-go thereby alleviating the low energy condition. The selection of a required 
OTT is a function of the amount of total energy at the start of the TAEM. 

The development of the TAEM trajectory for this report will cover only the nominal 
overhead entry point. The algorithms for the other approaches would be an extension of 
the basic nominal overhead algorithm, but adjusted based on the available energy. As a 
result, the trajectory equations will concentrate on the heading alignment cone 
trajectory profile. 
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Figure 5-3 HAC Configuration 






Figure 5-4 Optional TAEM HAC Trajectories 
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Figure 5-5 shows a graphical representation of the HAC trajectory parameters and 
variables. It is assumed the following initial HAC conditions and final conditions are 
known: 


- VELW1 = 

- HW1 = 

- VELNEP= 

- HNEP= 

- ^HAC= 

- <>HAC = 


Initial HAC velocity at WP1 (ft/sec) 
Initial HAC altitude (ft) 

Final HAC velocity at NEP (ft/sec) 
Final HAC altitude at NEP(ft) 
Constant HAC flight path angle 
Constant HAC bank angle 


Given the desired initial and final HAC conditions, the final fixed radius RF and initial 
radius at W1 can be calculated as follows: 

( 1 ) RF=(VHNEP**2)/(g*tan(<l> HAC )) where: VHNEP= horizontal velocity at NEP 


(2) RW1 =RF * (VHW1/VHNEP)**2 VHW1= horizontal velocity at W1 

The basic polar equation describing the spiral turn radius Is : 

RT=turn radius, for given 'P 

(3) RT= RF + R1* V P + R2*4'**2 where: RF= final turn radius ( V F=0) 

4' =HAC heading angle (deg) 
R1, R2= spiral coefficients 

From equation (3), the spiral coefficients R1, and R2 need to be determined to define the 
spiral turn radius RT. The primary purpose of the spiral is to enable a constant bank 
angle to be maintained during the HAC turn while modulating the velocity as given in 
equation (4). 

( 4 ) Vfioriz=SQRT[ 9 tan(<j>HAc)l 
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Figure 5-5 Definition of HAC Spiral Profile 




If the velocity were held constant during the HAC maneuver, the radius would be held 
constant generating a circular profile, resulting in a cylindrical shaped trajectory. 
However, PLS requires that the velocity decrease during the HAC maneuver while 
holding the bank angle constant; this requires a changing turn radius (RT). As the 
velocity decreases, the turn radius must also decrease resulting in a spiral profile. The 
spiral coefficients are used to shape the spiral to a desired profile as a function of the 
required velocity. 

The velocity constraints imposed by equations (1) and (2) may be used to solve for the 
spiral coefficients. In addition, there are altitude constraints imposed on the spiral 
trajectory. This constraint is imposed as result of the desired range during the HAC 
maneuver, as a function of flight path angle Yhac : 

( 5 ) RGhAC= (HW1-HNEP)/tan(Y HAC ) where: RG HAC =desired HAC range 

The resulting spiral range equation (RG'=predicted range) can be obtained by integrating 
the spiral equation (3) from 0 to the HAC heading angle 'f'(deg) : 

(6) RG'= 0 f RT (d^ *DTR) = [RF + RfY /2+ (R2*4'*‘2)/3]*4' *DTR 

where DTR = 1/57.3 (RAD/DEG) 

By setting the desired range RG HAC equal to the predicted range RG\ a solution for the 
required initial heading angle may be obtained. However, there is one equation with 

three unknowns. Experimentation showed that setting R1=0 and only using R2 yielded 
satisfactory trajectories. By setting R1=0, there are 2 unknowns, R2 and T 0 , and one 

equation. A quick iteration yields a solution determining both the spiral coefficient R2 
and the required initial heading angle 'f'o that will yield a spiral profile within given 

altitude and velocity constraints. Section 6.1 gives a solution to the required initial 
heading angle H'q and the spiral coefficient R2 for a sample TAEM trajectory. 

Having established the spiral coefficients, the desired trajectory profile command can be 
generated by equation (3), and the velocity command is defined by equation (4). 
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5.1.2 TAEM Altitude Control Module 


The TAEM altitude control module is identical to the A/L altitude control module. The 
altitude gains were developed for the A/L phase as a function of dynamic pressure. The 
TAEM dynamic pressure range is within the A/L range, therefore the gain schedules are 
sufficient for both phases. 


5.2 TAEM Lateral Guidance and Control System 

Lateral guidance is used during the TAEM heading alignment phase to track a nominal HAC 
(Heading Alignment Cone) spiral ground track profile. Ideally, if the vehicle maintains 
a constant bank angle in the absence of winds or external disturbance, the desired spiral 
ground track should be followed with minimum error. However, an external disturbance 
or deviation from the nominal velocity profile results in a change to the turn radius, and 
therefore ground track error. The lateral guidance system generates a roll command to 
eliminate the cross track error. The cross track error is the vehicle lateral AY error. 
The longitudinal AX error is controlled by the altitude and speed controller, and is 
covered in section 5.1.2 and 5.3. 

Our approach in the lateral guidance is to use the same cross track error and cross track 
rate law used in the A/L phase. The lateral guidance system used during the A/L phase as 
reported in section 3.3 consists of controlling the lateral displacement from the runway 
centerline, and the lateral velocity normal to the fixed centerline. The A/L control laws 
are easily adapted to the TAEM phase by using cross track error in place of the lateral 
displacement error and cross track rate error in place of the lateral velocity error. A 
block diagram of the lateral controller is given in figure 5-6. The bank angle generator 
is used to generate the proper HAC bank angle command as a function of the turn radius 
and the vehicle velocity. This relationship is defined by: 


(7) <f>HAC= tan-1 (V horiz 2/(g‘R T ) 




Figure 5-6 TAEM Lateral Control Law Block Diagram 


A simplified block diagram of the method used to calculated the cross track errors is 
given in figure 5-7. The detailed definition of each block shown in figure 5-7 is given 
in appendix H, which contains the full set of algorithms. A brief description of the cross 
track error and cross track rate error is presented below. 


X Y 



Figure 5-7 Block Diagram of Lateral Controller Cross 

Track Error Generator 
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A diagram showing the variables used to calculate the cross track error terms Is given In 
figure 5-8. The commanded ground track spiral profile and the actual ground track 
pro He is shown. The vehicle desired position is defined by the turn radius (RT) the 
commanded HAG ground track angle (Tc) defined in equation (3). and the commanded 
velocity (Vc) defined by equation (4). The resulting commanded radius vector is Rc . 

The actual vehicle position is defined by Ihe radius vector R. and the velocity vector V 

The commanded ground track angle (vc) is adjusted by the adaptive loop such thal for a 
corrected vc, 


( 8 ) dot product (Re • Vc)=0; 

that is, vc is adjusted to maintain Re perpendicular to Vc. Re is defined by: 

(9) Re=(X-Xc)i + (Y-Yc)j ; 
and the commanded velocity is defined by: 

(10) Vc= Xc i + Yc j . 


Substituting the above equations into the dot product equation (8), 


(11) dot product (Re • Vc)=0 = Xc(X-Xc) + Yc(Y-Yc) 

Equation (11) is used in the adaptive loop (see appendix H) to 
signal. As ¥e approaches zero, M's'Fc 

The cross track error becomes: 


create the 'Fe error 


(12) YCTerr=| R e |*sign(R-Rc)=SQRT[ (X-Xc)2 + <Y-Yc)2 ]*sign(R-R C ) ; 

and the cross track rate error becomes: 

(13) YCTerr=|(V x Vc)|/|Vc| = (X*Yc - Y*Xc)/SQRT(X C 2 + Yc 2) 
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The gains for the lateral steering were developed through nonlinear simulation analysis. 
Figure 5-9 shows the cross track error and bank angle command time histories as the 
vehicle captures the HAC spiral profile. The response is shown for combinations of 
cross track error gains and bank command gains. Nominal gain values selected were 
Ky=0.2 and K<t>=0.5. Using constant values for these gains, the PLS vehicle follows the 
desired HAC spiral profile within 50 feet throughout the heading alignment phase as 

shown in figure 5-10. The nominal error at the end of the phase is 33 feet in the X 
direction and 26 feet in the Y direction. 


85 



H 1 1 1 1 1 1 1 r 

i i i i i i i i i 

i 1 1 1 1 1 1 1 F 

i i i i i t i • 1 

i i i i i i\ i i i 

H f 1 1 1 hV- H 1 F 

I I I I I I \ < * 1 


+ - + - + - + 
flit 

+ - + - + - + 

(III 

+ - + - + - + 

lift 

+ - + - + - + ' 

I I 1 I 

+ ~ + ~ + ~ + 

till 

4 1 1 1“ 

ttti 

+ - + - + - + 

tiii 

+ - + - + - + 


4 b ir + 

i i i 

4 b ; ; + 

i i i 


+ - + 
t i 

+ - + 
I I 


I I I I I 

4 1 1 1 b 

i 1 t 1 1 1 

+ - + - 1 + - + - + 
I *\A T » t , 1 

- + -+ 

i -4. i^y i t\i i r-i 

|¥+r|f - + -?+ 


-+¥+ 


I I I I I I 

+ - + - + - + - + “ + " 

I I I I I I 

+ - + - + - + - + “ + - 
I I f t t t 

+ - + - + - + " + - + - 
1 I I > I I 

4 1 1 1 1 1 

( i t i i i 

_l 1 1 1 1 1 

i t i i i i 

+-+-+-+-+-+- 
i i i i t i 

+-+-+-+-+-+- 
i i i i i i 

4 1 1 1 1 1 

! i > I I i 

4 1 1 1 1 b - 

iiiiii/ 

4 1 1 1 b p4- 

\ i i t 1 o 1 

i i i i i ^ i 

4 1 1 1 b _ |4^> 

i i i i i 


t ^ c 

t i ^ 

— I b 

i i 

_ + __i 

i i 

- + - + 

I i 

“ 4 1 c 

i i * 

- 4 b 

I I 

— I 1 

I I 

-4-- + 

»» 1 1 

-&- + 
j. 

-/+ * 


(03SU) JJ913A 


-i r - “T t 1 

- i 1 r 

\ i i 

i i i i i 

H 1 1 1 1 

-*4 1 b ■ - 

i i i 

1 1 t 1 1 

4 1 1 1 1 

-4 b - 4- - - 

i i 

till! 

H 1 1 1 1 


1 1 i • 1 

4 1 1 1 ! 

4 -tit- 

1 1 i 1 1 


4 1 1 i 1 

1 1 i i 

1 1 1 1 1 
+ - + - + - + -H 

■ till 

1 1 . i i 

1 I i i 1 

4 1 1 1 1 

O 

- + + 

I 1 ! 1 1 

•jj ^ ! 

4 1 -1 1 1 

t | 1 

\ \ i i • 

4 1 1 1 1 

till 

| t i i 

1 1 1 

\ i i i * 

4 1* b ‘ 

\ i i ) 

i i i * i 

_i — i — i — i i - 

\+ - + - + ■ 
\ t i 

i i i i t ' 

i — i — i — i i 



4 b - 

4 b ' 

I I 

4 b 

f i 

4 b 

i i 

4 b 

i i 

4 b 

i i 

4 — b 

i ! 

4“ - 4- 

i i 

4 b 

t i 

4 b 

I 1 

H- - + 

I 1 

-h - + 

I t 


I 1 1 I t 

H 1 1 1 f- 

i i i i i 

-| 1 1 1 F • 

H 1 1 A ~ + 

i • i i\t 

4 1 1 hr + 

I 1 I 1 \ 1 

+-+-+-++ 

! I I I I I 

4 1 1 b * 4- 

i i i t i 

4 b -_4 1 b 

I ) O I i t 

■4 b ->4 b 4 4- 

i i it 

■4 — (■ Ah h - , b 

I 1 \ 1 1 1 

* 4 - H rV 4 — 14- 

i i 1 v i 1 

- 4- - 4- - + - V-H- 

- 4 1 1 ^ ~Ty 


DO W- 

(6ap) 3 <b 


Figure 5-9 TAEIi Lateral Control Performance For Different 

Lateral Gains 




a 

CO 



>- 

4-> 

a> 

Q 


/ifl+300 £ 


00+300 0 

OJ) X Eiiaa 


hO+300 <?- 


87 


Figures 10 TAEM Performance with Lateral Controller On Versus Off 



5.3 TAEM Speed Control 


The purpose of the PLS speed controller during the TAEM guidance phase is to control the 
vehicle to follow a predetermined speed profile. Speedbrake deflection is used to control 
the vehicle velocity. 

The A/L speed controller was adapted for use during the TAEM phase. The control law is 
identical {see figure 3-14), using an integral gain (Ksi) to reduce the regulator 
(velocity) steady state error, and gain Ks to reduce the velocity error. The speedbrake 

limit during TAEM is increased to 60 e . 

The velocity command (Vcmd) is based on the desired turn radius (equation 3) and the 
steady state bank angle command profile (<1>HAC> equ ation 7), resulting in equation (13): 


(13) V cm d= SQRT( Rt um * g * tan (4»HAC>] 

The 6- DOF simulation was used to evaluate the speed controller gains. The study showed 
that the gains used for the A/L phase result in sufficient speed control performance, and 
will therefore be used for the TAEM phase. 
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6.0 TAEM Six DOF Simulation Analysis 


In this section, the 6 DOF simulation is used in the mission design and evaluation of the 
TAEM performance. A TAEM parameter sensitivity study was conducted to define the 

nominal mission. The resulting mission design was verified through wind turbulence 
simulation. 


6.1 TAEM Parameter Sensitivity Study 

The 6-DOF simulation was used to evaluate the sensitivity of mission input parameters 
on the heading alignment phase final conditions. The results of these were used to select 
the nominal mission. The mission constraint is that an overhead approach to the runway 
be used to allow the pilot to view the runway before the landing. As a result, the heading 
angle changes during the heading alignment phase must be larger than 180 ? . 

The mission parameters analyzed, includes: 

• Y = HAC initial glideslope angle 

• <(> = steady state bank angle command 

• VELW1= velocity at way point 1 (start of HAC) 

• VELNEP= velocity at nominal entry point (end of HAC) 

• AH = altitude change during HAC 

Figure 6.1 shows heading angle change (A'F) sensitivity to the mission parameters. The 
sensitivity were run around the base point shown. The domain parameter is the steady 
state bank angle command with insufficient effect due to initial glideslope angle. 
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As a result of this study, the following mission input parameters will define the default 
TAEM mission: 

• Y = - 2 2 8 

•<{> = - 35 ® 

• VELW1= 900 ft/sec 

• VELNEP= 700 ft/sec 

• A H = 30000 ft (from 45000 to 1 5000 ft) 

With the above mission parameters, the HAC change in heading and spiral coefficient is- 

• A'F^ 187®. 

• R2=0.3505 


The guidance and control performance for the nominal mission are shown in figures 6-2 
through 6-5. Figure 6-2 contains 2 plots. Figure 6-2 shows the HAC commanded 
ground track (AY versus AX) and the resulting ground track profile; time histories for 

the ground track X and Y error are also shown. The position error and the end of the HAC 
flight is 33 feet in X and 26 feet in Y. 


The cross track error and rate error time histories for the lateral guidance system is 
plotted in figure 6-3 . Figure 6-4 shows the roll angle that results from the lateral 
error. The simulation starts with straight and level flight and proceeds to bank at the 
start of the HAC maneuver. As a result, there is an initial bank angle ramp to 
approximately 33 degrees. The vehicle bank angle response to the commanded bank angle 
shows well damped roll control performance. 

The flight path angle and angle of attack time histories are shown in figure 6-4 with 
altitude control time histories in figure 6-5 . The altitude error remains small and the 
error is only a couple feet at the end of the HAC phase. The altitude error is defined as 
the difference between the commanded and the actual altitude. Note, the commanded 
altitude at the end of the HAC flight is approximately 14620 feet, which is less than the 
targeted 15000 feet. This is due to the algorithm used to terminate the HAC flight phase. 
The algorithm uses the heading alignment angle as a test to terminate the HAC simulation; 

when the angle reaches zero, the simulation is terminated. This occurs at an altitude of 
14620 feet. 
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Figure 6-2 TAEM Simulation Ground Track Performance Results 
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Figure 6-3 TAEM Simulation Cross Track Error Results 
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Figure 6-5 TAEM Simulation Altitude Performance Results 
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6.2 


TAEM Performance With Wind Disturbance 


Simulation runs were also conducted using the wind luibulence model used for the A/L 
analysis (see section 4.3). The turbulence Id gaussian distribution and integral scale 

length is set according to the data in Appendix C. at an altitude of 9000 meters. These 

parameters are set at: 

• 0 y = 0£ = 17.3 ft/sec 

♦ L X = Ly= L Z= 1753 ft. 

Time history results for the HAC mission with wind turbulence, are presented in 
Appendix F A comparison of the final conditions with and without winds are shown in 
table 6-1. The wind turbulence has a very small effect on the final conditions. 


tafm End Condition 

Nominal 

Wind Turbulence 

Description 

AX err 

33 ft 

66 ft 

along track error 

AY err 

-26 ft 

-25 ft 

cross track error 

u 

14620 ft 

14612 ft 

altitude 


= 0 . 2 

- 0 . 2 

heading angle 


Table 6-1 TAEM End Conditions 
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6.3 TAEM Speed Control Performance 


Speed controller results from a nominal TAEM mission without wind is given in figure 
6-6. The initial velocity error oscillation is due to the initial speedbrake setting being 
30®. Despite this, the speed controller performed well up to 75 seconds. At this point, 
the vehicle velocity falls below the commanded velocity, and the speedbrake becomes 
fully retracted. 

The speed control performance can be easily improved through the re sizing of the turn 
radius constants used in the velocity command profile and bank command profile. This 
would have the effect of selecting an HAC spiral flight more consistent with the PLS drag 
characteristics. The speedbrake would operate between the full deflection and full 
retracted condition. 
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Figure 6-6 TAEM Simulation Speed Controller Results 
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7.0 Conclusion 


The NASA Langley 10 man PLS configuration is a very flyable configuration for 
performing a complete autoland mission similar to the previously studied 8 man CERV 
configuration. Acceptable landing performance can be achieved. If a manual backup or 
complete manual flight is required, the configuration is marginal and a higher L/D 
configuration is desired to ease the landing operation. The configuration has an L/D of 
3.2 and should be increased for manual flight. 

The autoland system performance was demonstrated using a non-linear 6-DOF 
simulation. The feasibility of using a double flare autoland profile was demonstrated 
with improved landing characteristics compared to the parabolic flare. The touchdown 
accuracy on sink rate and speed were within requirements. Adequate stability and 
control was demonstrated for wind turbulence and wind shear. 

Initial studies have been conducted on a portion of the TAEM system, including the 
heading alignment phase and prefinal phase. The feasibility of using a Shuttle type TAEM 
system for PLS was verified using a 6-DOF simulation; acceptable performance was 
demonstrated. Further studies will be required to complete the TAEM design including 
the design of the S-turn phase and acquisition phase. The purpose of these phases will be 
to transit from midcourse position and velocity uncertainties to the desired approach and 
landing initial conditions. 
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8.0 Recommendations 


Design studies have been conducted for the PLS Autoland system and simulation results 
demonstrate acceptable performance. To increase confidence in the approach, a number 
of configuration iterations and tailoring, along with control refinements and extension of 
the design could be pursued. Following are our suggestions for resolving some of the 
issues. 

The present PLS configuration has sufficient L/D for a pure autoland system, but for 
manual backup or manual flight, the L/D needs to be increased. One configuration 

suggestion would be to reduce the wing dihedral and increase the wing area to give 
additional lift required in the pitch plane. This would also reduce the C|p, which could 

help the trim during the decrab maneuver. However, because configuration change could 
adversely effect the high Mach number trim characteristic and yaw/roll stability, both 
high and low Mach regions would need to be studied. 

Landing in a cross wind has always been a problem with autolanding a lifting body. 
Vehicles with large C|p are difficult to trim at large sideslip angles. The decrab 

autopilot design in the CERV study (reference 7) was a complex design using an estimate 
of the trim state as feedforward command. Additional autopilot studies for the design of 
the decrab autopilot are needed. Relaxation of the decrab performance to only 3 fi of the 
runway centerline can be considered in the re-design. 

Areas where additional wind tunnel data is needed include: 

1) body flap deflection from zero to 40® to cover their use as speed brake and 

roll control; 

2) control surface hinge moment data is needed for actuator sizing; 

3) landing gear and ground effect for the final flare. 

Initial studies have been conducted on the last two segments of the TAEM system. An 
extension of this study is recommended to include the S-turn phase and acquisition 
phases of the TAEM system. A feasible study is needed to assure the PLS configuration 
can eliminate the re-entry guidance errors. This may be critical since the flight angle- 
of-attack capability of PLS is very low in this Mach region. 
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Also suggested is a feasibility study for the re-entry flight phase (440K ft to 60K ft 
altitude) to assure the stability and control characteristics, and to demonstrate re-entry 
phase accuracy. Initial 6-DOF simulations studies should be conducted using an 
integrated guidance and control design. 
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Appendix A 


AERODYNAMIC MODEL 


An aerodynamic model for the PLS vehicle was generated from the wind tunnel data 
provided. This data was obtained from NASA Langley Research Center. Force and moment 
equations of the model as a function of alpha and beta is given in figure A-1. Note, that 
this data is not a function of the control surface deflection. A simplification was used 
when generating this data, calculating the force and moment coefficient as a linear 
function of the control surface deflections. 

The nomenclature for the aerodynamic data is given in figure A-2; the definition of the 
control surface angles is given in figure 2-4. 

A trim program and derivative program were also developed. The derivatives are 
calculated around the trim conditions listed in section 3.1 . Sample derivatives at these 
trim conditions are given in figure A-3 . These derivatives were used to form the linear 
model for the autopilot point design analysis. 

The aerodynamic tabular data obtained from NASA is summarized as one large table . 
Figure A-4 shows the aerodynamic table format used to define the table. 
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AGCB PLS/CERV Subsonic Aerodynamic Model Description 


Build up .equations 


C x = Cx 0 (a.P) + AC X6c (a) 5e + ACx| 5 j(a) N + ACx 6f( (a) 5f + + ACx 5r (a) 5f- 
+ ACx| 5Af |(a) |SAf| + ACxj 5l |(a) N + 

C Y = C Yp P + ACy 5a (a) 5a + AC YgAf (a) 5Af + AC Ygr (a) 5r + 

Cz = Czo(a.P) + A Cz*( a > Be + AC 7 . 5r+ (a) 5C + ACz gr (a) 8f" 

Ci = C lp p + AC l5> (a) 5a + AC, 5Af (a) 5Af + AC, Sf (a) 8r + AC, p (a) 

+ ACj r (a) A + 

C m = C^a.P) + AC mj . c (a) 5e + AC m j 5a j(a) |5a| + AC m ^(a) 5f* + AC m5r (a) 8f 
+ AC^/a) |5n| + AC mq (a) ^ ° 

C n = C no (a,p) + AC n5a (a) 8a +AC n5Af (a) 5Af + AC n{ . f (a) 8r + AC„ p (a) 

+ AC nr (a) A + 
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AGCB PLS/CERV Subsonic Aerodynamic Model Description 



Note: a and are shown as representative independent 
variables; refer to the cross reference table above to 
determine actual independent variables. 


Figure A-4 Aerodynamic Table Format 
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CXO 


Columns 

1 thru 

6 




-I . 00000+01 

-3.8089D-02 

-4 . 2403D-02 

-4. 59110-02 

-4 . 9690D-02 

-5.3195D-02 

-6.5021D+00 

-5 . 00 61D-02 

-5.3242D-02 

-5 . 5828D-02 

-5.8872D-02 

-6.17 98D-02 

- 3 .03240+00 

-5. 98630-02 

-6.19210-02 

-6.3593D-02 

-6.5907D-02 

-6.8258D-02 

3 . 8120D-01 

-6.53 1 8D-02 

-6.6271D-02 

-6.7043D-02 

-6.06390-02 

-7 . 04 2 6D-02 

3 . 7 1 12D+00 

-6.5707D-02 

-6.5582D-02 

-6.54760-02 

- 6 . 637 2D-02 

-6. 7608D-02 

6. 9309D+00 

-6.14 23D-02 

-6 . 0256D-02 

-5. 9301D-02 

-5 . 9S20D-02 

-6.0223D-02 

1 .0014D+01 

-5.3571D-02 

-5.1406D-02 

-4 . 9 638D-02 

-4 . 9208D-02 

-4 . 9401D-02 

1 . 2 936D+0 1 

-4.35800-02 

-4 . 04 69D-02 

- 3 . 7 932D-02 

-3.68870-02 

-3.6597D-02 

1. 56740+01 

-3 . 2884D-02 

-2.80070-02 

-2 . 5628D-02 

-2.4008D-02 

-2 . 32 64D-02 

1.82050+01 

-2 . 2690D-02 

-1. 70730-02 

-1 . 3947D-02 

-1.1795D-02 

-1 . 0632D-02 

2.05080+01 

-1 . 3850D-02 

-8.2879D-03 

-3 . 7544D-03 

-1.1177D-03 

4 .26000-04 

2. 25660+01 

-6.8265D-03 

-5 . 9827D-04 

4 . 4778D-03 

7.5473D-03 

9.4315D-03 

2.4362D+01 

-1 . 72 40D-03 

5.08550-03 

1 .06350-02 

1 . 4082D-02 

1 . 62 63D-02 

2.5881D+01 

1 . 6300D-03 

8.9312D-03 

1 .48810-02 

1 .8640D-O2 

2. 10810-02 

2 . 7 112D+01 

3 . 5833D-03 

1 . 1283D-02 

1.7557D-02 

2. 1582D-02 

2.421 9D-02 

2 . 8043D+01 

4 . 5551D-03 

1 . 255 6D-02 

1. 90760-02 

2 . 3297D-02 

2.6088D-02 

2 . 8 668D+01 

4 . 9439D-03 

1. 31470-02 

1 . 9832D-02 

2.41850-02 

2.7079D-02 

2 . 8982D+01 

5.0557D-03 

1 .3361D-02 

2 . 0128D-02 

2.45470-02 

2 . 74 93D-02 

2 . 9061D+01 

5.0748D-03 

1 .34050-02 

2 . 01 93D-02 

2 . 4628D-02 

2 . 7587D-02 

2.9375D+01 

5.11480-03 

1 .3547D-02 

2.0418D-02 

2 . 4919D-02 

2 . 7 930D-02 

3 . 0000D+01 

, 5.01800-03 

1 . 3 653D-02 

2 .06890-02 

2.5321D-02 

2 . 84 36D-02 

O.OOOOD+OO 

p-1 . OOOOD+Ol 

-7 . 5001D+00 

-5.4643D+00 

-3 . 8402D+00 

-2.5622D+00 

Columns 

7 thru 12 




-5.6026D-02 

-5.8150D-02 

-5 . 9701D-02 

-5 . 8 150D-02 

-5. 6026D-02 

-5 . 3195D-02 

-6.4 172D-02 

-6 . 5923D-02 

-6 . 7 143D-02 

-6.5923D-02 

-6. 41720-02 

-6. 1798D-02 

-7.01800-02 

-7 . 1561D-02 

-7 . 2 453D-02 

-7 . 1561D-02 

-7.0180D-02 

-6 . 8250D-O2 

-7 . 1901D-02 

-7 . 2 9 18D-02 

-7.34000-02 

-7 . 2918D-02 

-7 . 1901D-02 

-7. 04260-02 

- 6 . 8648D-02 

-6 . 93 10D-02 

-6.95650-02 

—6 . 9310D-02 

-6.8648D-02 

-6.7608D-02 

-6 . 0843D-02 

-6.1161D-02 

-6 . 1 112D-02 

-6 . 11 61D-02 

-6 . 0843D-02 

-6 . 0223D-02 

-4 . 9619D-02 

-4 . 9608D-02 

-4 . 92670-02 

-4 . 9608D-02 

-4.96190-02 

-4 . 9401D-02 

-3 . 6433D-02 

-3. 61100-02 

-3.5493D-02 

-3.6110D-02 

-3 . 6433D-02 

-3 . 6597D-02 

-2.2743D-02 

-2.21280-02 

-2. 12530-02 

-2 . 2 1 28D-02 

-2.2743D-02 

-2 . 3264D-02 

-9.7806D-03 

-8 . 8964D-03 

-7.7812D-03 

-8 * 8964D-03 

-9.7806D-03 

-1 . 0632D-02 

1.5786D-03 

2.7086D-03 

4 . 04 15D-03 

2 . 7086D-03 

' 1 . 5786D-03 

4 .26000-04 

1 . 0853D-02 

1.2202D-02 

1 .3730D-02 

1 .2202D-02 

1 . 0853D-02 

9.431 5D-03 
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1. 79190-02 

1.9460D-02 

2 . 1158D-02 

1 . 9460D-02 

1.7919D-02 

1 . 62 630-02 

2 . 2 935D-02 

| 2 . 4 638D-02 

2 . 64790-02 

2 .46380-02 

2 . 2935D-02 

2 . 1081D-02 

2 . 6234D-02 - 

' 2 . 80 6BD-02 

3 . 0025D-02 

2.8068D-02 

2. 62340-02 

2 . 42 19D-02 

2 . 8224D-02 

3.01580-02 

3.22030-02 

3 . 0158D-02 

2 . 8224D-02 

2 . 6088D-02 

2 . 92 97D-02 

3.1297D-02 

3 . 3401D-02 

3. 12970-02 

2.9297D-02 

2.7079D-02 

2. 97520-02 

1 3 . 1786D-02 

3.3919D-02 

3.1786D-02 

2 . 9752D-02 

2 . 74 93D-02 

2 . 9857D-02 

1 3.18 99D-02 

3 . 4040D-02 

3 . 18 99D-02 

2.98570-02 

2.7587D-02 

3 . 0240D-02 

3.231 6D-02 

3 . 4487D-02 

3. 23160-02 

3. 02400-02 

2 . 7 930D-02 

3 . 0828D-02 

3.2970D-02 

3 . 5200D-02 

3.29700-02 

3 . 0820D-O2 

2 . 84 36D-02 

-1 . 5538D+00 

"-7 . 300 ID-01 

0.0000D+00 

77300 10-0 1 

”1 .55300+00 

2 ."56220 + 00 

7 

8 

7 

to 

ft 



Columns 

13 thru 

16 


-4 . 9690D-02 

-4 . 591 ID-02 

-4 . 2 403D-02 

-3 . 8089D-02 

-5 . 8872D-02 

-5.5828D-02 

-5.32420-02 

-5 . 006 ID-02 

-6 . 5907D-02 

-6.3593D-02 

-6.192 ID-02 

-5.98630-02 

-6 . 8639D-02 

-6 . 7043D-02 

-6 . 6271D-02 

-6.53 18D-02 

-6 . 6372D-02 

-6.54 7 6D-02 

-6.5582D-02 

-6.5707D-02 

-5 . 9520D-02 

-5. 9301D-02 

-6. 02560-02 

-6. 1423D-02 

-4 . 9208D-02 

- 4 . 9638D-02 

-5 . 1406D-02 

-5.357 ID-02 

-3.6887D-02 

-3 . 7 932D-02 

-4 . 04 69D-02 

-4 . 3580D-02 

-2 . 4OO0D-O2 

-2 . 5 620D-O2 

-2.88870-02 

-3 . 2884D-02 

-1. 17950-02 

-1.394 7D-02 

-1 .7873D-02 

-2 . 2 690D-02 

-1. 1177D-03 

-3 . 7544D-03 

-8.2879D-03 

-1 . 3850D-02 

7 . 5473D-03 

4 . 4778D-03 

-5 . 9827D-04 

-6 . 82 65D-03 

1. 40820-02 

1 .06350-02 

5 . 0855D-03 

-1 .7240D-03 

1.8648D-02 

1 . 4881D-02 

8 . 93 12D-03 

1 . 6300D-03 

2 . 1582D-02 

1 .7557D-02 

1 . 1283D-02 

3 . 5833D-03 

2.3297D-02 

1.9076D-02 

1.2556D-02 

4 . 555 ID-03 

2 . 4 185D-02 

1 . 9832D-02 

1. 31470-02 

4 . 9439D-03 

2 . 4 547D-02 

2 . 0 1 28D-02 

1. 33610-02 

5 . 0557D-03 

2 . 4 628D-02 

2.01 93D-02 

1.3405D-02 

5.0748D-03 

2.49190-02 

2.04 18D-02 

1 . 3547D-02 

5 . 1148D-03 

2 . 5321D-02 

2. 06890-02 

1.3 6530-02... 

5 . 0188D-03 

3 . 8402D+00 

5.4643D+00 

7 . 5001D+00 

1 .00000+01 
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czo 


Columns 
-1 .00000+01 
-6. 5021D+00 
-3 .03240+00 
3 . 8 120D-01 
3.7112D+00 
6. 9309D+00 
i.-, 1.0014D+Q1 
1. 29360+01 
1 . 56740+01 
1 . B205D+01 
2 .05080+01 
2 . 2566D+01 
2 .43620+01 
2 . 5881D+01 
2 .7112D+01 
2. 80430+01 
2 . 8 668D+01 
2 . 8 982D+01 
2 . 9061D+01 
2 . 9375D+0 1 
3 . 0000D+01 
0 . 0000D+00 

Columns 
5. 1316D-01 
3 . 501 4D-01 
2.0603D-01 
6 . 9389D-02 
-6. 4611D-02 
-1.9613D-01 
-3 . 2272D-01 
-4.4 10 ID-01 
“5 . 4795D-0 1 
-6. 4147D-01 
-7 .20740-01 
-7 .8597D-01 


i thru 
4 . 4846D-01 
2.9980D-01 
1. 69940-01 
4 .73140-02 
-7.3011D-02 
-1 . 9131D-01 
-3.0524D-01 
-4.1153D-01 
-5 . 0723D-0 1 
“5 . 90360-01 
-6. 6017D-01 
-7 . 1 695D-0 1 
-7.6179D-01 
-7. 96150-01 
-8.2156D-01 
-8.3938D-01 
-8 . 5067D-01 
-8. 56140-01 
-8.5749D-01 
-8.6279D-01 
-8.72 99D-01 
-1 .OOOOD+Ol 


6 

4 . 6077D-0 1 
3 .07850-01 
1.7378D-01 
4.7002D-02 
-7.7371D-02 
-1 . 9959D-01 
-3.1726D-01 
-4.27100-01 
-5 . 2 613D-0 1 
-6.12340-01 
-6.8494D-01 
-7 . 4 423D-01 
-7 . 9125D-0 1 
-8.27460-01 
-8 . 5436D-01 
-8.7331D-01 
-8.85360-01 
-8 . 9121D-01 
-8 . 9265D-01 
-8 . 9834D-0 1 
-9 . 0930D-0 1 
-7 . 5001D+00 


4 .7619D-01 
3 . 1981D-01 
1 . 8230D-01 
5 . 2 148D-02 
-7 . 5520D-02 
-2 . 0092D-01 
-3 . 2165D-01 
-4 . 3438D-01 
-5 . 3 612D-01 
-6.2483D-01 
-6 . 9971D-01 
-7.6104D-01 
-8 . 0984D-01 
-8.4755D-01 
-8. 75670-01 
-8 . 9554D-01 
-9.08210-01 
-9. 1437D-01 
-9.15890-01 
-9 . 2189D-01 
-9.3346D-01 
-5 . 4 643D+00 


5. 20510-01 
3 . 5608D-01 
2 . 1059D-01 
7 . 2578D-02 
-6.. 27550-02 
-1.95570-01 
-3.2339D-01 
-4.42850-01 
-5 . 5088D-01 
-6. 4542D-01 
-7 . 2561D-01 
-7 . 91670-01 


4 . 9109D-01 
3.31950-01 
1 . 9170D-01 
5.8852D-02 
-7.14460-02 
-1. 99390-01 
-3.22550-01 
-4.37590-01 
-5 . 4 149D-01 
-6 . 3220D-01 
-7 . 0890D-01 
-7 . 7185D-01 
-8.2207D-01 
-8. 6098D-01 
-8 . 9007D-01 
-9. 1068D-01 
-9.23840-01 
-9. J025D-01 
-9.31830-01 
-9. 3808D-01 
-9. 5015D-01 
-3 . 8402D+00 


5 . 1316D-Q1 
3 . 5014D-01 
2.06030-01 
6.93890-02 
-6 . 4611D-02 
-1 . 9613D-01 
-3 . 2272D-01 
-4 , 4 101D-01 
-5 . 4795D-01 
-6, 4147D-01 
-7 . 2074D-01 
-7.8597D-01 


5 . 0353D-01 
3.4222D-01 
1 . 99820-01 
6.4844D-02 
-6.7523D-02 
-1 . 97 47D-01 
-3 . 2254D-01 
-4 . 3940D-01 
-5. 4500D-01 
-6.37280-01 
-7.1541D-01 
-7 .7964D-01 
-8 . 3098D-01 
-8.7083D-01 
-9 . 00 69D-01 
-9 . 2187D-01 
-9 3542D-01 
-9.42030-01 
-9. 43660-01 
-9. 5010D-01 
-9. 6256D-01 
-2.5622D+00 


5 . 0353D-01 
3.4222D-01 
1.9982D-01 
6.48440-02 
-6.75230-02 
-1.9747D-01 
-3 . 2254D-01 
-4 . 3940D-0 1 
-5.4500D-01 
-6 . 3728D-01 
-7. 15410-01 
-7. 79640-01 


7 thru 12 
5 . 2051D-01 
3.56080-01 
2.10590-01 
7.2578D-02 
-6.2755D-02 
-1 . 9557D-0 1 
-3.23390-01 
-4 .42850-01 
-5.5088D-01 
-6 . 4 542D-01 
-7 . 25 61D-01 
-7 . 9167D-01 


J . £. OJHU-U1 

3.60670-01 
2 . 139SD-01 
7 . 4724D-02 
-6.17910-02' 

r3.2466D-0h 
-<T4tl?i)-d‘l 
-5. 54170-01 
-6 . 4 960D-01 
-7 . 3061D-01 
-7.9740D-01 


AGCB PLS/CERV Subsonic Aerodynamic Model Function Tables February 15, 1990 


-8 . 3819D-01 

-8 . 4460D-01 

-8.50970-01 

-0 . 4460D-01 

-8.38190-01 

-8 . 3098D-01 

-8.7079D-O1 

-8.85810-01 

-8 . 927 ID-01 

-8 . 8581D-01 

-8.78790-01 

-8.7083D-01 

-9 . 0925D-01 

-9.16760-01 

-9.24 10D-01 

-9 . 1676D-01 

-9.09250-01 

-9 . 0069D-01 

-9.30090-01 

-9.38780-01 

-9 . 4645D-01 

-9.30780-01 

-9. 30890-01 

-9.21870-01 

-9.4475D-01 

-9.52890-01 

-9.6078D-01 

-9 . 5289D-01 

-9.44750-01 

-9 . 3542D-01 

-9.5 1 5QD-01 

-9 . 5977D-01 

-9 . 6777D-01 

-9.5977D-01 

-9 . 5150D-01 

-9 . 4203D-01 

-9 . 53 18D-01 

-9. 61 47D-01 

Sj £9510-01 

-9 . 6147D-01 

-9.5318D-01 

-9.4366D-01 

-9 . 5977D-01 

-9 . 68190-01 

-9. 7634*5-01 

-9 . 6819D-01 

-9 . 5977D-01 

-9 . 5010D-G1 

-9.7253D-01 

-9 . 8121D-01 

-9 . 8957D-Q1 

-9 . 8121D-01 

-9.7253D-01 

-9. 62560-01 

-1. 55380+00 

-7 . 3001D-01 

0.00000+00 

7 . 3001D-01 

1 . 5538D+00 

2 . 5 622D+00 

Columns 

13 thru 16 




4 . 9109D-01 

4.76190-01 

4 . 6077D-01 

4 . 404 6D-01 



3 . 3 1 95D-01 

3 . 1981D-01 

3. 07850-01 

2 . 9980D-01 



1 . 9 170D-01 

1 . 82 30D-01 

1 .73780-01 

1 . 6994D-01 



5 . 8852D-02 

5.2148D-02 

4.70020-02 

4 . 7314D-02 



-7 . 1446D-02 

-7.5520D-02 

-7 . 7371D-02 

-7 .3011D-02 



-1. 99390-01 

-2.0092D-01 

-1 . 9959D-01 

-1.91310-01 



-3 . 2255D-01 

-3.21650-01 

-3.1726D-01 

-3 . 0524D-01 



-4 . 37 59D-01 

-4 . 3438D-01 

-4.2710D-01 

-4 . 1153D-01 



-5. 4149D-01 

-5. 3612D-01 

-5 . 2 613D-01 

-5. 07230-01 



-6 . 3220D-01 

-6 . 2483D-01 

-6 . 1234D-01 

-5 . 9036D-01 



-7 . 08 90D-01 

-6. 9971D-01 

-6. 84940-01 

-6.60170-01 



-7.7185D-01 

-7 . 6104D-01 

-7. 44230-01 

-7. 16950-01 



-8 . 2207D-0 1 

-8 .09840-01 

-7. 91250-01 

-7 . 617 9D-01 



-8.6098D-01 

-8.4755D-01 

-8 . 2746D-01 

-7. 96150-01 



-8 . 9007D-01 

-8. 7567D-01 

-8.5436D-01 

-0 . 2156D-01 



-9. 1068D-01 

-8 . 9554D-01 

-0.73310-01 

-0 . 3938D-0 1 



-9.23840-01 

-9.0821D-01 

-8. 85360-01 

-8 . 5067D-Q1 



-9 . 3025D-01 

-9. 1437D-01 

-8 . 912 ID-01 

-8. 56140-01 



- 9 . 3183D-01 

-9. 1589D-01 

-8 . 9265D-01 

-8.5749D-01 



-9 . 3808D-01 

-9.21890-01 

-8 . 9834D-01 

-8 . 6279D-01 



-9 . 5015D-01 

-9.3346D-01 

-9.0930D-01 

-8 . 7299D-Q1 



3 . 8402D+Q0 

5 . 4 643D+00 

7 .50010+00 

1 .OOOOD+Ol 
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CMO 

Columns 
-1 .OOOOD+Ol 
-6.50210+00 
-3 . 0324D+00 
3 . 8 120D-0 1 
3 . 7112D+00 
6. 9309D+00 
1 . 0014D+01 
1.2936D+01 
1 . 5674D+01 
1 -02O5D+O1 
2 . 0508D+01 
2 . 2566D+01 
2 . 43 62D+01 
2 . 5801D+O1 
2 . 7 112D+0 1 
2. 8043D+01 
2 . 8668D+01 
2 . 0982D+O1 
2 . 9061D+01 
2. 9375D+01 
3 .OOOOD+Ol 
0 . OOOOD+OO 

Columns 
4 . 74280-02 
3 . 8 96 ID-02 
3. 18950-02 
2 . 4787D-02 
1 .7495D-02 
1 .0422D-02 
4 . 02 74D-03 
-1.3872 D-03 
-5 . 7006D-03 
-8 . 9001D-03 
-1 . 1006D-02 
-1 .20440-02 


1 thru 
3 . 8066D-02 
3 . 1390D-02 
2 . 6100D-02 
2 . 0739D-02 
1.5152D-02 
9.72670-03 
4 . 9106D-03 
9.9185D-04 
-1 . 9203D-03 
-3 . 8242D-03 
-4 . 7507D-03 
-4 . 7 3 50D-03 
-3.84 42D-03 
-2 . 2300D-03 
-1.674 6D-04 
1 . 9542D-03 
3 . 6979D-03 
4 . 6813D-03 
4 . 9396D-03 
6 .0210D-03 
8 . 4281D-03 
-1 .OOOOD+Ol 


6 

4 .0272D-02 
3.3066D-02 
2 . 7250D-02 
2 . 1372D-02 
1 . 5280D-02 
9.36700-03 
4 .0845D-03 
-2 . 7702D-04 
-3 . 6039D-03 
-5 . 8 913D-03 
-7. 1 668D-03 
-7.46300-03 
-6. 8442D-03 
“5.4 602D-03 
-3.58410-03 
-1 . 6035D-03 
4 . 5395D-05 
9. 8128D-04 
1 .2276D-03 
2 . 2 615D-03 
4.5739D-03 
-7.50010+00 


4 .2887D-02 
3 , 5250D-02 
2 . 9006D-02 
2.2706D-02 
1 - 6204D-02 
9 . 8937D-03 
4 .2300D-03 
-4 . 9204D-04 
-4 . 1567D-03 
-6 .7563D-03 
-8.3161D-03 
-8 . 86620-03 
-8 . 4690D-03 
-7 . 2724D-03 
-5 . 5481D-03 
-3 . 6825D-03 
-2 . 1107D-03 
-1 . 2 135D-03 
-9.7689D-04 
1 .82290-05 
2 . 2535D-03 
-5 . 4 643D + 00 


4 . 8290D-02 
3 . 9649D-02 
3 . 2409D-02 
2.51300-02 
1 - 7672D-02 
1 . 0439D-02 
3 . 8902D-03 
-1 . 6703D-03 
-6 . 1204D-03 
-9.44 62D-03 
-1 . 1667D-02 
-1 .2808D-02 


4 . 4914D-02 
3. 6932D-02 
3.0346D-02 
2.37110-02 
1. 68810-02 
1 .0254D-02 
4 . 2865D-03 
-7.23130-04 
-4.6572D-03 
-7.5060D-03 
-9.2925D-03 
-1 . 0045D-02 
-9.8248D-03 
-8.7777D-03 
-7. 1745D-03 
-5 . 4006D-03 
-3.8903D-03 
-3.024 ID-03 
-2.79520-03 
-1 . 8309D-03 
3.4283D-04 
-3.84020+00 


4 . 7428D-02 
3.89610-02 
3. 1895D-02 
2.4767D-02 
1 .7495D-02 
1 .0422D-02 
4 .0274D-03 
-1.3872D-03 
-5.7006D-03 
-8.90010-03 
-1 . 1006D-02 
-1 . 204 4D-02 


7 thru 12 
4 . 82 90D-02 
3 . 9649D-02 
3 . 2409D-02 
2 . 5130D-02 
1.76720-02 
1 .0439D-02 
3 . 8 902D-03 
-1.6703D-03 
-6 . 1204D-03 
-9.4462D-03 
-1.1667D-02 
-1 . 2808D-02 


4 . 9085D-02 
4 .0290D-02 
3 . 2897D-02 
2.54 67D-Q2 
1 . 7861D-02 
1 . 04 85D-02 
3 . 8003D-03 
-1.8895D-03 
-6 . 4 607D-03 
-9. 8985D-03 
- 1 . 2 221 D -02 
-1 .3453D-02 
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-1.2072D-02 
-1 . 1236D-02 
-9 .8029D-03 
-8 . 1580D-03 
-6.734 4D-03 
-5.911 6D-03 
-5. 6936D-03 
-4 .7729D-03 
-2 . 6857D-03 
-1 . 5538D+00 

Columns 
4 . 4914D-02 
3 . 6932D-02 
3 . 034 6D-02 
2.3711D-02 
1 . 688 ID-02 
1 . 0254D-02 
4 . 28 65D-03 
-7 . 2313D-04 
-4.6572D-03 
-7 . 50600-03 
-9.29250-03 
-1 .0045D-02 
-9.8248D-03 
-8.77770-03 
-7 . 1745D-03 
-5.40060-03 
-3 . 8903D-0 3 
-3.024 ID-03 
-2. 7952D-03 
-1 .8309D-03 
3.4283D-04 
3 . 8402D+00 


-1 .29260-02 
-1.21 65D-02 
-1 . 0794D-02 
-9 , 1 953D-03 
-7 . 8029D-03 
-6. 99580-03 
-6.7817D-03 
-5.8767D-03 
-3 . 8208D-03 
-7.3001D-01 


-1 .3650D-02 
-1 .2957D-02 
“1 . 1640D-02 
-1 .0083D-02 
-8.71820-03 
-7 . 9250D-03 
-7. 7144 D-03 
-6.82320-03 
-4 . 7 94 9D-03 
0. OOOOD+OO 


-1 .2926D-02 
-1 .2165D-02 
-1 . 07 94D-02 
-9. 1953D-03 
-7 . 8029D-03 
-6 . 9958D-03 
-6.7817D-03 
-5 . 8767D-03 
-3 . 8208D-03 
7 . 30010-01 


3 . 8066D-02 
3 . 1 390D-02 
2 . 6100D-02 
2.073 9D-02 
1 . 5152D-02 
9.72670-03 
4 . 9106D-03 
9 . 9185D-04 
-1 . 92030-03 
-3 . 82 4 2D-03 
-4 . 75070-03 
“4 . 7350D-03 
-3.84 4 2D-03 
-2.23000-03 
-1 . 6746D-04 
1 . 95420-03 
3 . 6979D-03 
4 . 6813D-03 
4 . 9396D-03 
6 . 0210D-03 
8 .42810-03 
1 . OOOOD + Ol 


13 thru 
4 . 28870-02 
3 . 5250D-02 
2 . 9006D-02 
2 . 2 706D-02 
1 . 6204D-02 
9 . 8937D-03 
4 .2300D-03 
-4 . 9204D-04 
-4 . 1567D-03 
-6 . 7563D-03 
-8.3161D-03 
-8.86620-03 
-8 . 4690D-03 
-7 . 2724D-03 
-5. 5481D-03 
-3. 6825D-03 
-2 . 11070-03 
- 1 . 2 135D-03 
-9 . 7689D-04 
1 .8229D-05 
2.25350-03 
5 . 4 643D + 00 


16 

4 . 0272D-Q2 
3.3066D-02 
2.7250D-02 
2 . 1372D-02 
1 .5280D-02 
9.36780-03 
4 .08450-03 
-2 . 7702D-04 
-3 . 6039D-03 
-5 .891 3D-0 3 
-7 . 1668D-03 
-7 . 4630D-03 
-6 . 8442D-03 
-5 . 4 602D-03 
-3 . 58410-03 
-1 - 6035D-03 
4 . 5395D-05 
9 . 8128D-04 
1.22760-03 
2.26150-03 
4.573 9D-03 
7 . 50010+00 


-1.2072D-02 - 

-1 . 12360-02 - 

-9 . 8029D-03 -j 
-8 . 1580D-03 
-6 . 7344D-03 -J 
-5.9116D-03 
-5 . 6936D-03 -4 

-4 . 7729D-03 -j 
-2.68570-03 -1 

1.5538D+00 2 
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4 . 6361D-02 
3 . 8 108D-02 
3. 1254D-02 
2 . 4354D-02 
1.72660-02 
1 . 038 9D-02 
4 . 1833D-03 
-1.05270-03 
-5.198 9D-03 
-8 .24360-03 
-1 .02090-02 
- 1 . 1121 D -02 
-1 . 10390-02 
- 1 . 01100-02 
-8 . 6021D-03 
-6. 9003D-03 
-5.43850-0 3 
-4 . 5965D-03 
-4 .3737D-03 
-3 . 4338D-03 
-1 . 3084D-03 
2 . 5 622D + 00 


4 . 63610-02 
3 . 8108D-02 
3 . 1254D-02 
2 . 4354D-02 
1.72660-02 
1 . 0389D-Q2 
4 . 1833D-03 
1 .0527D-03 
5 . 1989D-03 
8 . 243 6D-03 
1 . 0209D-02 
1 . 1121D-02 


February 15, 1990 


1 . 1039D-02 
lUiOl 10D-02 
1 . 602 ID-03 
i. 90030-03 
- . 4385D-03 
. 5965D-03 
.37370-03 
. 4338D-03 
. 3084D-03 
. 5622D+00 


Table (continued) 



AGCB PLS/CERV Subsonic Aerodynamic Model Function Tables 


February 15, 1990 


CYB 

- 1 .24210-02 
CLNO 


CLLB 

-7 .891 9D-03 


Columns 

1 thru 

6 




-1 .00000+01 

-2.3000D-02 

- 1 . 1500D-02 

-4 .6000D-03 

0. 00000 + 00 

4 . 6000D-03 

-5 . 0000D+00 

-2 . 8000D-02 

-1 .4000D-02 

-5. 60000-03 

0. OOOOD+OO 

5 . 6000D-03 

0 . OOOOD+OO 

-4. OOOOD-02 

-2 . 0000D-02 

-1 .0000D-02 

0. OOOOD+OO 

1. 0000D-02 

5 . OOOOD+OO 

-3. OOOOD-02 

-1 . 4000D-02 

-5 . 000QD-03 

0. OOOOD+OO 

5. 00000-03 

1 .OOOOD+Ol 

-3. OOOOD-02 

-1 .8000D-02 

-8 . OOOOD-03 

0. OOOOD+OO 

8. 0000D-03 

1 .5000D+01 

-3 . 0000D-02 

-1 . 5000D-02 

-6.0000D-03 

0. OOOOD+OO 

6. 0000D-03 

2 .OOOOD+Ol 

-1.2000D-02 

-1 . 0000D-02 

-7 .00000-03 

0 .OOOOD+OO 

7. 0000D-03 

2 . 5000D+01 

-9. 9000D-03 

-5 . 3000D-03 

-1 .7000D-03 

0. OOOOD+OO 

1 . 7000D-03 

3. OOOOD+Ol 

-2 . 0000D-04 

2 . 4000D-03 

2 . 9000D-03 

0. OOOOD+OO 

-2 . 9000D-03 

0. OOOOD+OO 

-1. OOOOD+Ol 

-5. OOOOD+OO 

-2 . OOOOD+OO 

0. 00000+00 

2. OOOOD+OO 

Columns 

7 thru 

8 




1 . 1500D-02 

2 . 3000D-02 





1 . 4000D-02 

2 . 8000D-02 





2 . OOOOD-02 

4 . 0000D-02 





1 . 4000D-02 

3. 00000-02 





1.8000D-02 

3 .0000D-02 





1.5000D-02 

3. 00000-02 





1 .00000-02 

1 . 2000D-02 





5 . 3000D-03 

9. 90000-03 





-2 . 40000-03 

2 .00000-04 





5. 00000+00 

1 .OOOOD+Ol 
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CXDE 


CZDE 


-1. OOOOD+Ol 

4 .4375D-04 

-1. OOOOD+Ol 

-4 . 1348D-03 

-6.2900D+00 

2 . 6802D-04 

-6 . 2900D+00 

-4 . 6625D-03 

-3 . 1200D+00 

2 . 0531D-04 

-3. 1200D+00 

-4 . 9652D-03 

-4 , 3000D-01 

1 .86770-04 

-4.3000D-01 

-5 . 1227D-03 

1 .8400D+00 

1 .81770-04 

1 . 8400D+00 

-5 . 1868D-03 

3 . 7500D+00 

1.78 17D-04 

3 . 7500D+00 

-5 . 1918D-03 

5 . 3600D+00 

1.721 6D-04 

5 . 3600D+00 

-5 . 1605D-03 

6 . 7300D+00 

1. 6325D-04 

6 . 7300D+00 

-5.1074D-03 

7.9200D+00 

1 .5182D-04 

7 . 9200D+00 

-5 . 0410D-03 

0 . 9900D+00 

1,381 6D-04 

8 . 9900D+00 

-4 . 9644D-03 

1 .OOOOD+Ol 

1 . 2206D-04 

1. OOOOD+Ol 

-4.8770D-03 

1 . 1010D+01 

1 . 02 66D-04 

1 . 10100+01 

-4 . 7743D-Q3 

1 .2080D+01 

7.840 ID-05 

1 . 2080D+01 

-4.648 ID-03 

1.3270D+01 

4^69flBD-Qi 

1.3270D+01 

-4 . 485 9D-03 

1 . 4640D+01 

/T78688D-06> 

1 .4640D+01 

-4 . 2690D-03 

1 . 6250D+01 

-T.T989D-05 

1 . 6250D+01 

-3 . 9700D-03 

1 . 8160D+01 

-1 . 2859D-04 

V. 8160D+01 

-3.5487D-03 

2 . 0430D+01 

-2 . 2926D-04 

2 . 0430D+01 

-2.9448D-03 

2 . 3120D+01 

-3.5455D-04 

2.3120D+01 

-2 . 0667D-03 

2 . 6290D+01 

-4 . 9384D-04 

2 . 6290D+01 

-7.7336D-04 

3 .OOOOD+Ol 

-6 . 1043D-04 

3 .OOOOD+Ol 

1 . 1547D-03 
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CMDE 

- 1 . 0000D+01 
-6 . 2900D+00 
-3 . 1200D+00 
-4 . 3000D-01 
1 .8400D+00 
3 . 7500D+00 
5 . 3600D+00 
6. 73000 + 00 
7 . 9200D + 00 
8 . 9900D + 00 
1 .OOOOD+Ol 
1 . 1010D+0 1 
1 . 2080D+01 
1.3270D+01 
1 . 4 640D+01 
1 . 6250D+01 
1 .81600+01 
2 . 0430D+0 1 
2 . 3120D+01 
2 . 62 90D+01 
3 .OOOOD+Ol 


-1.54 96D-05 
-1 .'7144D-03 
- 1 . 82 92D-03 
-1 .8953D-03 
-1 . 9228D-03 
-1 . 92 33D-03 
-1 . 9065D-03 
-1 .8793D-03 
-1 . 84 60D-03 
-1 .8083D-03 
- 1 . 7 660O-03 
-1 .7172D-03 
-1 . 6585D-03 
-1.5852D-03 
-1 . 4905D-03 
-1 . 3 659D-03 
-1 . 20140-03 
-9. 85380-04 
-7 .07 6 6D-04 
-3 . 6398D-04 
3 . 3662D-05 
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CXDDE 

-1.00000+01 -9.03770-04 

-6 . 2900D+00 -9 . 0377D-04 

-3.12000+00 -9.0377D-Q4 

-4.3000D-01 -9.0377D-04 

1 . 8400D+00 -9 . Q377D-04 

3.75000+00 -7.9201D-04 

5.36000+00 -6 . 3599D-04 

6 . 7300D+00 -5.96040-04 

7 . 9200D+00 -5. 22710-04 

8 . 9900D+00 -4.6317D-04 

1. OOOOD+Ol -4 . 1415D-04 
1.1010D+01 -3.7324D-04 

1 . 2080D+01 -3 . 3965D-04 

1.3270D+01 -3.147 ID-04 

1.4640D+01 -3 . 0240D-04 

1.6250D+01 -3 . 0932D-04 

1 . 8160D+01 -3.43270-04 

2 . 04 30D + 0 1 -4.0784D-04 

2 . 312GD+0 1 -4 . 8833D-04 

2 . 6290D+0 1 -5 . 1985D-04 

3. OOOOD+Ol -3.2150D-04 


Table A-1 PLS Aerodynamic 


CMDDE 

-1. 0000D+01 2 . 07 97D-04 

-6 . 2900D+00 2 . O797D-04 

-3. 12000+00 2.0797D-04 

-4 . 3000D-01 2.0797D-04 

1.84000+00 2.0797D-04 

3.7500D+00 2 . 2992D-04 

5.3600D + 00 2 . 277 2D-04 

6 . 7300D+00 2.1762D-04 

7 . 9200D+00 2 . 0628D-04 

8 . 9900D+00 1.9612D-04 

1. 00000+01 1.8790D-04 

1. 10100+01 1 . 8 199D-04 

1.20800+01 1 . 7 910D-04 

1 . 3270D+01 1.8070D-04 

1 . 4 640D+01 1 . 893 6D-04 

1 . 6250D + 01 2 . 088 ID-04 

1 .8160D+01 2 . 4305D-04 

2 . 04 30D+01 2 . 9302D-04 

2 . 3120D+01 3.47650-04 

2.6290D+01 3.63310-04 

3. OOOOD+Ol 2 . 2 101D-04 
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CYDDE 

-1 . 0000D+0 1 

3 . 3156D-03 

CLLDDE 

-6.2900D+00 

3 . 3156D-03 

-1 .0000D+01 

-3 . 1200D+00 

3. 3156D-03 

-6 . 2900D+00 

-4 .30000-01 

3 . 3156D-03 

-3 . 1200D+00 

1 .8400D+00 

3.3156D-03 

-4 . 3000D-01 

3 . 7500D+00 

3 .2514D-03 

1 .8400D+00 

5.36000+00 

3 . 1819D-03 

3.75000+00 

6 . 7300D+00 

3 . 11280-03 

5 .36000+00 

7 . 9200D+00 

3 . 0458D-03 

6.73000+00 

8. 99000+00 

2 . 9806D-03 

7. 9200D+00 

1 .00000+01 

2 . 9148D-03 

8 . 9900D+00 

1. 1010D+01 

2 . 8454D-03 

1 .0000D+01 

1 .2080D+01 

2 . 76810-03 

1 .1010D+01 

1 . 3270D+01 

2 . 6779D-03 

1.20800+01 

1 . 4640D+01 

2 . 5691D-03 

1.3270D+01 

1 . 6250D+01 

2 . 4352D-03 

1 .46400+01 

1.8160D+01 

2.26890-03 

1 . 6250D+01 

2 .0430D+01 

2 . 0627D-03 

1.81600+01 

2.31200+01 

I . 8090D-03 

2 . 0430D+01 

2 . 62 90D+01 

1 . 5010D-03 

2 . 3 120D+01 

3 . 00000+01 

1.13360-03 

2 . 62900 + 01 


3 . 0000D+01 


February 15 


2 . 5618D-03 
2 . 5618D-03 
2 . 5618D-03 
2.561 00-03 
2 . 56 18D-03 
2.561 ID-03 
2.541 6D-03 
2.51190-03 
2 .4767D-03 
2.43800-03 
2 . 3955D-0J 
2.3473D-03 
2.29020-03 
2 . 21 90D-O3 
2.12990-03 
2 . 012 9D-03 
1.85870-03 
1 - 6552D-03 
1 .3874D-03 
1.0378D-03 
5.8611D-04 
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CLNODC 

-1 .00000+01 
-6.29000+00 
-3 . 1200 D +00 
-4 . 3000D-01 
1 . 84000+00 
3 . 75000+00 
5.36000+00 
6.73000+00 
7.92000+00 
8 . 9900D+00 
1 . 00000+01 
1 . 1010 D +0 1 
1 .2Q80D+01 
1 . 327QD+01 
1.46400+01 
1 . 6250D+01 
1 .8160D+01 
2 . 0430D+01 
2 . 3 120D+0 1 
2 . 6290D+01 
3 . 0000D+0 1 


-2.817BD-03 
-2.8178D-03 
-2 . 0170D-O3 
-2 , 8178D-03 
-2.81700-03 
-2.81120-03 
-2 . 7696D-03 
-2.71390-03 
-2.65350-03 
-2 .59180-03 
-2.52840-03 
-2 . 4 61 3D-03 
-2.3872D-03 
-2 . 3026D-03 
-2 . 203 70-03 
-2 . 08760-03 
-1 . 9524D-03 
-1. 79730-03 
-1.6213D-03 
-1.4173D-03 
-1. 1579D-03 


, 1990 


1990 
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AGCfl PLS/CKRV Subsonic 


Aerodynamic Model 


funci. ion Tables 


CXDBFP 


February 15, 1990 


-1 . OOQOD+0 1 

-6 . 2900D+00 
-3. 1200D+00 
“4 . 30000-01 
1 .8400D+00 
3 . 7500D+00 
5 . 36000+00 
6. 7300D+00 
7 . 92000+00 
8 . 9900D+OQ 
1 .OOOOD+01 
1 . 10100+01 
1 .2080D+01 
1 . 3270D+01 
1 • 4640D+01 
1.62500+01 

1 . 8160D+01 
2 . 0430D+01 
2 . 3 120D+01 
2 . 62900+01 
3 . OOOOD+Ol 


-2 . 7210D-04 
'1 - 9054D-04 
~ 1 • 4 61 3D-04 
-1.31330-04 
-1 ■ 36320-04 
-1 . 5274D-04 
-1.74730-04 
-1 . 98750-04 
-2.2313D-04 
-2.47440-04 
-2. 7210D-04 
-2 . 98070-04 
-3.26550-04 
-3 . 5878D-04 
-3 . 9559D-04 
-4 .3673D-04 
-4 • 7941D-04 
-5.15570-04 
-5.2697D-04 
-4 .76100-04 
-2 . 9061D-04 


CZDBFP 


-1 .0000D+01 
-6 . 2900D+00 
-3 . 1200D+00 
-4 . 3000D-01 
1 . 8400D+00 
3 . 75000+00 
5 . 3600D+00 
6 . 7300D + 00 
7 . 92000+00 
8 . 9900D+00 
1 .00000+01 
1 . 1010D+01 
1 .2080D+01 
1.3270D+01 
1.4640D+01 
1 • 6250D+01 
1 .8160D+01 
2 . 04 30D+01 
2 .3120D+01 
2. 6290D+01 
3.00000+01 


-3. 7100D-03 
-3 .75120-03 
-3. 77200-03 
-3.77870-03 
-3 . 77 64D-03 
-3.7689D-03 
-3 . 7587D-03 
-3.74730-03 
-3.7353D-03 
-3.72300-03 
-3 . 7 100D-03 
-3.6958D-03 
-3 . 6794D-03 
“3.65960-03 
-3. 6348D-03 
-3. 6033D-03 
-3 . 56260-03 
-3. 5103D-03 
-3 . 4437D-03 
-3 . 3 607D-03 
-3.2607D-03 


AGCB PLS/CERV Subsonic 


Aerodynamic Model 


Function 


Tables 


CMDBFP 


-1 . 00000 + 01 
-6 . 2900D+00 
-3. 1200D+00 
-4 , 30Q0D-Q1 
1 .8400D + 00 
3 . 7500D+00 
5 . 3 600D + 00 
6. 7300D+00 
7 . 9200D+00 
8 . 9900D+00 
1 . OOOOD+Ol 
1- 10100+01 
1 .2080D+01 
1 .32700+01 
1 . 4 640D+01 
1 • 62500+01 
1 • 016OD+O1 
2 . 0430D + 01 
2 . 3 120D+01 
2 . 6290D + 0 1 
3 . OOOOD+Ol 


-9. 8909D-04 
-9.8953D-04 
-9.8958D-04 
-9.89570-04 
-9.8957D-04 
-9 . 8958D-04 
-9 . 895 6D-04 
-9. 8951D-04 
-9.8 94 ID-04 
-9. 8928D-04 
-9.89090-04 
-9 . 8802D-O4 
“9.884 4D-04 
-9.8788D-04 
-9. 8699D-04 
-9 . 0557D-O4 
-9 . 8320D-04 
-9.7 92 ID-04 
-9.7236D-04 
-9 . 6049D-04 
-9.3986D-04 


February 15, 1990 
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CXDBFN 

-1. OOOOD+Ol 
-6.2 900D+00 
-3 . 1200D+00 
-4 . 3000D-01 
1 .840004-00 
3 . 7500D4-00 
5 . 3600D4-00 
6.73000+00 
7 . 9200D+00 
8.99000+00 
1 .0000D+01 
1.10100+01 
1 . 2080D+01 
1.32700+01 
1 .46400+01 
1 . 6250D+01 
1.81 60D+01 
2 . 0430D+01 
2 . 3120D+01 
2 . 62 90D+01 
3 . 0000D+01 


4 . 22 41D-04 
5.521 ID-04 
4 .75400-04 
3 . 4820D-04 
2. 31840-04 
1.44020-04 
8 . 4 652D-05 
4 . 8043D-05 
2 .80930-05 
2 .01840-05 
2 . 1651D-05 
3.1762D-05 
5 . 1630D-05 
8. 41160-05 
1.3359D-04 
2 . 0504D-04 
3.01460-04 
4 . 1752D-04 
5 .25440-04 
5 . 4 634D-04 
2 . 9495D-04 


CZDBFN 

-1.0000D+01 
-6 . 2900D+00 
-3. 12000+00 
-4 .30000-01 
1 . 8400D+00 
3.75000+00 
5 . 3600D+00 
6 . 7300D+00 
7 . 9200D+00 
8.99000+00 
1 .OOOOD+Ol 
1.10100+01 
1 .2080D+01 
1.32700+01 
1 . 4640D+01 
1 . 6250D+01 
1 .81600+01 
2.04300+01 
2 . 3120D+01 
2 . 6290D+01 
3 .OOOOD+Ol 


-3.44610-03 

-3.7456D-03 

-3 . 7862D-03 

-3.7261D-03 

-3 . 6387D-03 

-3.5547D-03 

-3 . 484 ID-03 

-3 . 4282D-03 

-3.38480-03 

-3 . 3510D-03 

-3.3243D-03 

-3.3033D-03 

-3.2876D-03 

-3 . 2785D-03 

-3.2789D-03 

-3.2939D-03 

-3.3305D-03 

-3 . 3952D-03 

-3.48930-03 

-3.5947D-03 

-3 . 647 ID-03 


AGCB PLS/CERV Subsonic Aerodynamic Model Function Tables 


February 15, 1990 


CMDBFN 

- 1 . 00000+01 
-6.2900D+00 
-3 . 1200D+00 
-4 . 3000D-01 
1 .8400D+00 
3. 75000+00 
5.36000+00 
6.73000+00 
7 . 9200D+00 
8. 9900D+00 
1. 0000D+01 
1.1010D+01 
1.20800+01 
1 .32700+01 
1 .4640D+01 
1 . 6250D+01 
1 .8160D+01 
2 . 0430D+01 
2 . 3120D+01 
2 . 62 90D+01 
3. OOOOD+Ol 


-1 . 06I2D-03 
-1 . 1362D-03 
-1 .1276D-03 
-1 . 0931D-03 
-1 .0568D-03 
-1 . 0271D-03 
-1 .0054D-03 
-9. 9077D-04 
-9.8142D-04 
-9.7 603D-04 
-9 . 7369D-04 
-9.74 10D-04 
-9.77510-04 
-9 . 84Q0D-04 
-9 . 9742D-04 
-1 .0172D-03 
-1 .0459D-03 
-1 . 0830D-03 
-1.1222D-03 
-1. 1427D-03 
-1.0927D-03 
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CXDDBF 


CYDDBF 


-1 . 0000D+0 1 
-6.2 900D+00 
-3. 1200D+00 
-4 .3000D-01 
1 .8400D+00 
3 . 7500D+00 
5 . 3 600D+00 
6 . 7300D+00 
7 . 9200D+00 
8 . 9900D+00 
1 . OOOOD+Ol 
1 .1010D+01 
1 .20800+01 
1 .3270D+01 
1 .4640D+01 
1. 6250D+01 
1 .81600+01 
2 . 0430D+01 
2 . 3120D+01 
2 . 6290D+01 
3 .0000D+01 


-2.7645D-04 
-4 . 4161D-04 
-5. 3697D-04 
-5.9619D-04 
-6 . 3 61 ID-04 
-6. 6506D-04 
-6.8730D-04 
-7 . 05 15D-04 
-7.200 6D-04 
-7 . 3307D-04 
-7.4502D-04 
-7 . 5 665D-04 
-7 . 6853D-04 
-7.011 ID-04 
-7 . 9449D-04 
-8 . 0814D-04 
-8 .2024D-04 
-8 . 2 632D-04 
-8 . 16560-04 
-7 . 70700-04 
-6 . 4848D-04 


-1 .OOOOD+Ol 
-6 . 2900D+00 
-3 . 1200D+00 
-4 .30000-01 
1 .8400D+00 
3.7500D+00 
5 . 3 600D+00 
6 .73000+00 
7. 92000+00 
8 .99000+00 
1.00000+01 
1 . 1010D+01 
1.20800+01 
1.32700+01 
1 .4640D+01 
1 . 6250D+01 
1 . 0 160D+01 
2 . 0430D+01 
2 . 3120D+01 
2 . 6290D+01 
3 .OOOOD+Ol 


4.351 9D-04 
2 . 8 17 7D-04 
1 .5067D-04 
4 . 3358D-05 
-4 . 244 6D-05 
-1 . 1033D-04 
-1 . 6404D-04 
-2 . 0701D-04 
-2 . 42 17D-04 
-2.7204D-04 
-2 . 9066D-O4 
-3.2376D-04 
-3.4868D-04 
-3.7438D-04 
-4 . 0138D-04 
-4 . 2967D-04 
-4.5867D-04 
-4.8726D-04 
-5 . 14 17D-04 
-5.38910-04 
-5 . 6388D-04 
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CLLDDBF - CLNDDBF 


-1 .OOOOD+Ol 

7.072 ID-04 

-1 .OOOOD+Ol 

-1.46440-04 

-6. 29000+00 

8.252 ID-04 

-6 . 2900D+00 

-1 . 1726D-04 

-3.1200D+00 

7 . 9730D-04 

-3 . 1200D+00 

-8 . 5288D-05 

-4 . 3000D-01 

7 . 5307D-04 

-4 . 3000D-01 

-5.5830D-05 

1 .8400D+00 

7 . 1215D-04 

1 .04000+00 

-3 . Q833D-05 

3 . 7500D+00 

6.8035D-04 

3 . 7500D+00 

-1.05380-05 

5 . 3600D+00 

6.5775D-04 

5.36000+00 

5.5843D-06 

6 . 7300D+00 

6 . 42 60D-04 

6.73000+00 

1 . 8334D-05 

7 . 9200D+00 

6 . 32 94D-04 

7.92000+00 

2.8532D-05 

8.9900D+00 

6.2721D-04 

8 . 9900D+00 

3.6903D-05 

1. OOOOD+Ol 

6.244 6D-04 

1. OOOOD+Ol 

4 . 4040D-05 

1 . 1Q10D+01 

6.2425D-04 

1. 1010D+01 

5 . 0367D-05 

1 . 2080D+0 1 

6 . 2 67 3D-04 

1 . 2080D+01 

5 . 6124D-05 

1 . 327QD+0 1 

6. 3253 D- 04 

1 . 3270D+01 

6.13010-05 

1 . 4640D+01 

6.4269D-04 

1 .4640D+01 

6 . 5560D-05 

1 . 6250D+0 1 

6.5838D-04 

1 . 6250D+01 

6.8085D-05 

1 . 8 1 60D+01 

6 . 8003D-04 

1.81600+01 

6.7385D-05 

2 . 0430D+01 

7.0528D-04 

2 . 0430D+01 

6. 1023D-05 

2.3120D+01 

7.244 6D-04 

2.3120D+01 

4 .5296D-05 

2.6290D+01 

7 . 1 152D-04 

2 . 6290D+01 

1.4897D-05 

3. OOOOD+Ol 

6.064 ID-04 

3 .OOOOD+Ol 

-3. 7330D-05 
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CXDR 


CYDR 


-1 .0000D + 01 
-6.2900D+00 
-3 . 1200D+00 
“4 . 3000D-0 1 
1.8400D+00 
3 . 7500D+00 
5 . 3600D+00 
6 . 7 300D+00 
7. 9200D+00 
0. 9900D+00 
1 . OOOOD+Ol 
1. 1010D+01 
1 .2080D+01 
1 . 3270D+01 
1 . 4 640D+01 
1 . 6250D+01 
1.8160D+01 
2 . 0430D+01 
2 , 3120D+01 
2 . 6290D+01 
3. OOOOD+Ol 


-6 . 0020D-04 
-6 . 0020D-04 
-6 . 002QD-04 
“6 . 0020D-04 
-6 . 0020D-04 
-6 . 1220D-04 
-6. 2276D-04 
-6 . 3S28D-04 
-6. 5020D-04 
-6. 6744D-04 
-6.8723D-04 
-7.1042D-04 
-7. 3843D-04 
-7 . 7309D-04 
-8 . 1616D-04 
-8 . 6803D-04 
-9.2463D-04 
-9.7056D-04 
-9.6462D-04 
-0. 1120D-04 
-3 . 0594D-04 


-l .00000+01 
-6.2900D+00 
-3 . 1200D+00 
-4 . 3000D-01 
1 .84000+00 
3 . 7500D+00 
5 . 3600D+00 
6 . 7300D+00 
7 . 9200D+00 
8 . 9900D+00 
1 .0000D+01 
1 . 1010D+01 
1 .2080D+01 
1 .3270D+01 
1 .4640D+01 
1 . 6250D+01 
1 .8160D+01 
2.0430D+01 
2 . 3120D+01 
2 . 6290D+01 
3 .OOOOD+Ol 


1 . 8950D-03 
i .8950D-03 
1 . 8950D-03 
1 . 8950D-03 
1 . 8950D-03 
1 . 9730D-03 
2.0622D-03 
2 . 15 13D-03 
2 . 2363D-03 
2.3173D-03 
2 . 3964D-03 
2 . 4770D-03 
2 . 5628D-03 
2 . 657 ID-03 
2 . 7620D-03 
2 . 6766D-03 
2 . 994 7D-03 
3 . 1000D-03 
3 . 1591D-03 
3 . 1104D-03 
2.8473D-03 
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CLLDR 


-1 .OOOOD+Ol 

2. 1894D-04 

-6 . 2900D+Q0 

2 . 1894D-04 

-3.1200D+00 

2 . 18 94D-04 

-4 . 3000D-01 

2 . 18 94D-04 

1 . 8400D+Q0 

2 . 1894D-04 

3.7500D+00 

2 . 3902D-04 

5 . 3600D+00 

2 . 6954D-04 

6.7300D+00 

3.00820-04 

7 . 9200D+00 

3 . 2 961D-04 

8 . 9900D+00 

3.554 6D-04 

1 .OOOOD+Ol 

3.7899D-04 

1.1010D+01 

4 . 0 1 10D-04 

1 .2080D+01 

4 . 2251D-04 

1 .3270D+01 

4.4355D-04 

1 . 4640D+01 

4 . 6404D-04 

1 . 6250D+01 

4 . 8353D-04 

1 .01600+01 

5 . 0232D-04 

2 . 0430D+01 

5.2464D-04 

2 . 3120D+01 

5.6622D-04 

2 . 6290D+01 

6. 7108D-04 

3 .OOOOD+Ol 

9. 4552D-04 


CLNDR 

-1. OOOOD+Ol -1 .26840-03 
-6 . 2900D+00 -1.2604D-O3 

-3 . 1200D+00 -1 .2684D-03 

-4 . 3000D-01 -1.2684D-03 

1 . 8400D+00 -1.2684D-03 

3 . 7500D+00 -1.2832D-03 

5 .36000+00 -1.3093D-03 

6.7300D+00 -1 . 3377D-03 

7 . 9200D+00 -1 . 3652D-03 

8.9900D+00 -1.3908D-03 

1. 0000D+01 -1 . 4 150D-03 

1 . 1010D+01 -1 . 4384D-03 

1 . 2080D+0 1 - 1 . 4618D-Q3 

1.3270D+01 -1 . 4854D-03 

1 . 4 640D+01 -1.5086D-03 

1 . 6250D+01 -1.5295D-03 

1 . 8160D+0 1 *-1.54 4 6D-03 

2 . 0430D+0 1 -1 . 5487D-03 

2.3120D+01 -1 . 5368D-03 

2 . 6290D+01 -1.5087D-03 

3. OOOOD + Ol -1.4815D-03 


Table A-1 PLS Aerodynamic Model Function Table (continued) 



AGCB PLS/CERV Subsonic Aerodynamic Model Function Tables 


February 15, 1990 


CMQ 

0 . OOOOD+OO 
5 . OOOOD+OO 

1 . OOOOD+O 1 
1 . 2000D+0 1 
1 .4000D+01 
1. 6000D+01 
1 . 8000D+01 

2 .0000D+01 
2 . 2000D+01 
2 .50000+01 
3 . OOOOD+0 1 

CLP 

0 .0000D+00 
5 .OOOOD+OO 
1.0000D+01 
1 . 2000D+0 1 
1 .4000D+01 
1 . 6000D+01 
1 .8000D+01 
2 . 0000D+01 
2 . 2000D+01 
2.5000D+01 
3.0000D+01 


-2 . 03J6D-01 
-1. 58360-01 
-1.1635D-01 
-1.7 63 6D-01 
-1 . 7636D-01 
- 1 . 7 53 6D-01 
-1.74 3 6D-01 
-2 . 5237D-01 
-2 . 9937D-01 
-5 . 7 13 6D-01 
-8.Q036D-01 


-4 . 98 1 ID-01 
-6 . 001 ID-01 
-3 . 981 ID-01 
-5 . 7 91 ID-01 
-4.251 ID-01 
-4 . 9911D-01 
-6 . 4911D-01 
-6. 1911D-01 
-8.101 ID-01 
-8.721 ID-01 
-2 . 4 10 1D+00 


CLR 

0. OOOOD+OO 
5. OOOOD+OO 
1 . OOOOD+Ol 
1 . 2000D+01 
1 . 4000D+01 
1 . 6000D+01 
1 .8000D+01 
2 .OOOOD+Ol 
2 . 2000D+01 
2 . 5000D+01 
3 .OOOOD+Ol 


4 .9613D-01 
5 . 9013D-O1 
8 . 4013D-01 
8 . 0613D-01 
6.8613D-01 
5 . 7213D-01 
6 . 0313D-01 
8 . 9913D-01 
9.5713D-01 
1 .3571D+00 
4 . 4901D+00 
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CNP 


CNR 


0 . OOOOD + OO 
5. OOOOD+OO 
1. OOOOD+Ol 
1 . 2000D+01 
1 . 4000D+01 
1 . 6000D+01 
1 .8000D+01 
2 .OOOOD+Ol 
2 . 2000D+01 
2 .50000+01 
3. OOOOD+Ol 


3 . 80ST7D-01 
3.5297D-01 
2 . 1897D-01 
2 . 43 97D-01 
1 . 7897D-01 
1 . 6697D-Q1 
2 . 0097D-01 
2 . 2 197D-01 
3 . 0697D-01 
3.7897D-01 
8 . 4997D-01 


0. OOOOD+OO 
5. OOOOD+OO 
1 .OOOOD+Ol 
1 . 2000D+01 
1 . 4000D+01 
1 . 6000D+01 
1 .8000D+01 
2. OOOOD+Ol 
2 . 2000D+01 
2 . 5000D+01 
3 .OOOOD + Ol 


-7 .9434D-01 
-8 . 3734D-0 1 
-9 . 2134D-01 
-9.2234D-01 
-8 . 6734D-01 
-9 . 2234D-01 
-9.4634D-01 
-1.0073D+00 
-1 . 0903D + 00 
-1 .2863D + 00 
-2 . 04 13D + 00 
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Appendix B 


6D0F NON-LINEAR SIMULATION 


The computer resources used in the 6DOF simulation are shown in figure B-1 
The non-linear equations of motion are hosted on the AD100. The autopilot and other 
slower processes are hosted on the Harris 1000. The PLS simulation runs three times 
real time. Figure B-2 shows the PLS simulation functional block diagram. The major 
blocks consist of the autopilot, guidance. 6DOF equations of motion, actuator models rate 
gyro models, and the gain computer. The autopilot operates at a 20 ms cycle, with a 10 
ms computational delay. Autopilot gains are updated from the gain computer every 120 
ms. Perfect knowledge of a. 6, and 0 is assumed and used as feedback signals. Autopilot 
fm commands are fed to a mixer set equations to obtain the proper command to the 
actuators. The actual fin positions are- then un-mixed before going to the aerodynamic 
model. The actuators are modeled as a second order filter with a 10 hz bandwidth and 0.5 
damping, and a rate limit of 200 degrees/sec. 

The guidance operates al a 60 ms cycle. Perfect knowledge of the spacecraft 
position and inert, al velocities is assumed and used as feedback signals to the guidance, 
ongitudinal and lateral guidance designs are described in more detail in sections 3.1. 
Rate gyro sensors are modeled as second order filters with a bandwidth of 60 hz and 0 4 
damping. The turbulence wind model is described in more detail in Appendix C 
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Figure B-2 PLS Simulation Block Diagram 



Appendix C 


WIND TURBULENCE MODEL 


C.1 Introduction 


number of efforts have been directed toward generating mathematical models of wind 
urbulence near the surface of the earth. References 3-6 all recommend the use of the Von 

Herman rsotroprc power spectra over the Dryden spectra, apparently due to their greater 
accuracy at high frequency. 

The Von Harman spectra provide a complex spacial distribution of wind velocity vectors 
w ic must be converted to a temporal distribution, by aircraft speed and flight path for use in a 

Z TT B ThS * ,he airCra "’ in - - - -ocrlies. must a, so be 72 

she “ 3 P,DP0SeS 10 aPPro ’" ma,e ,he Von Kama " distributions by driving 

a thir^lT^ Tr- " P ' 0P0SeS 3 *« *>' iongitudina, direction and 

basT ,h 81 ^ d ' reC,i0nS ' Re,erCnCe 4 ate ° I"*— shaping titters 

based on the landrng conditions of a transport aircraft, forthis approximation. I, uses a first-order' 

wet for the longrtudinal direction and a second-order liner for the lateral and vertical directions 
This is the model that has been implemented in our simulation. 

C.2 Model Description 

Tbe form ot the wind turbulence model shaping tillers used in the flight simulation is taken 
tmm Reference 4. stock diagrams ol the continuous troms o. the filters, in terms ol displacement 
a ong the II, gh. path, are shown in figures C-t and C-2. In the simulation these filters are 
discretized and the aircraft speed is used to convert lime steps to distance steps The 

~r hTr de ' emi " ed ,mm ,he ai ' C,a " ”' i " 9 Span ,b) and,he ch ° se " W^tra, 

cate leng,h (L). whrch determines Ihe Irequency response of the liners. The chosen variance 

eve, (sigma) of .he wind componems. along wilh the Integra, scale length, determines .he level o, 
the input white noise to drive the gust spectra. 

The choice o, appropriate integral scale length and turbulence variance is dilticult since 
these vary with altrtude and measured data show large variations. Figures C-3 and C-4 trom 
e erence 4 provide an rndication ol the probabililies and measured variations in these values 

rgure C-5 tram Re, erence 6 gives altitude pretties o, these values as chosen tor use in Space 
Shuttle simulations. F 
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Figure C-1 Turbulence Model lor X-Axis 






Figure C-2 Turbulence Model for Y and Z-Axis 
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Rgure C-3 Probability Distribution of the RMS Turbulent Velocity 
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Figure C-4 Variation of the Turbulence Scale with Height 


VARIATION OF STANDARO DEVIATION 
AND LENGTH SCALE WITH ALTITUDE* 


ALTITUOE 

(m) 


10 

20 

30 

40 

50 

60 

70 

80 

90 

100 

200 

304.8 

400 

500 

600 

700 

762 

800 

900 

1524 

2000 

3048 

4000 

5000 

6096 

7000 

8000 

9144 

10000 

20000 


STANDARD DEVIATION 
OF TURBULENCE 


INTEGRAL SCALES 
OF TUR8ULENCE 


*1 (m/sec i a 2(m/ sec) 


3.72 3.53 

3.95/4.37 3.95/4.37 
4.39 4.39 

4.39 4.39 


4.39/5.704.39/5.70 
5.70 5.70 

5.70 5.70 

5.70/5.79 5.70/5.79 
5.79 5.79 

5.79/5.52 5.79/5.52 
5.52 5.52 

5.52 5.52 

5.52/5.27 5.52/5.27 
5.27 5.27 

5.27 5.27 

5.27/4.22 5.27/4.22 
4.22 4.22 

6.01 6.01 


a 3 (m/sec 


1.15 

1.46 

1.71 

1.89 

2.05 

2.19 

2.32 

2.43 

2.54 

2.64 

3.38 

3.95/4.39 

4.39 
4.39 
4.39 
4.39 

4.39/5.70 


LjOn) 

L 2 (m) 

L 3 ( m ) 

jn 

11 

5 

33 

19 

11 

43 

28 

17 

52 

35 

23 

61 

42 

29 

6R 

49 

35 

75 

56 

41 

82 

63 

47 

89 

69 

53 

95 

75 

59 

149 

134 

123 


196/300 

300 

300 

300 

300 

300/533 

533 

533 

533 

533 

533 

533 

533 

533 

533 

533 

533 

533 

6691 


190/300 

300 

300 

300 

300 

300/533 

533 

533 

533 

533 

533 

533 

533 ; 

533 

533 

533 

533 

533 

6691 


Double entries for a tabulated altitude indicate a step 
change in standard deviation or Integral scale at that altitude. 


Figure C-5 Variation of Standard Deviation and Length Scale With Altitude 
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Appendix D 


Simulation Results of Approach and Landing Nominal Mission 

Without Wind 


The nominal approach and landing mission was simulated without wind. The conditions of 
this simulation are defined in section 4.2 .The results are plotted in the following 
figures. These results show the nominal autoland performance . 
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Appendix E 


Simulation Results of Approach and Landing Nominal Mission 

With Wind Turbulence 


The nominal approach and landing mission was simulated with wind turbulence. The 
conditions of this simulation are defined in section 4.3 . The results are plotted in the 
following figures. These results show excellent touchdown performance. 
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Appendix F 


Simulation Results of a TAEM Nominal Mission 
With Wind Turbulence 


A nominal TAEM mission was simulated with wind turbulence. The conditions of this 
simulation are defined as follows. The turbulence la gaussian distribution and integral 

scale length is set according to the data in Appendix C. at an altitude of 9000 meters. 
These parameters are set at: 

•ox=a Y =az= 17.3 ft/sec 
* L X= l Y= l Z= 1753 ft. 


The results are* plotted in the following figures. 
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Appendix G 


Autoland Approach and Landing Trajectory 
Generation Algorithms 

The following section contains a block diagram of the trajectory generator used to 
generate guidance commands during the autoland A/L phase. These algorithms also define 
the second order exponential transition equations used in guiding the vehicle from the 

end of the pre-flare maneuver ( from altitude of H2C to H2), and during the final flare 
(altitude H3 to touchdown). 

The definition of the autoland A/L trajectory is covered in section 3.2.1 . 
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Initial value: 

ACIR=Pull-up g's = 0.5 


HSTAR-RN1 

H DOTC - H S T AR * S DOT 

HC-HC+ ( HDOTC+HDOTCP ) *TAU60/2 

HDOTCP-HDOTC 

PULUP-0 




Figure G-l PLS TRAJECTORY ALGORITHMS BLOCK DIAGRAM 
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Appendix H 


Cross Track Error/ Cross Track Rate Error 
Generation Algorithms 


The following section contains a block diagram (see figure H-1) of the cross track error 
and cross track rate error generator used to generate bank angle guidance commands 
during the TAEM phase. These algorithms serve to define the functions given in figure 
5-7. A unique feature is the adaptive loop which calculates d('Pe)dt. This error 
function serves to adapt the trajectory command such that when, d(*f'e)dt.=0 , then the 
commanded angular rate [d('Pc)dtJ equals the actual rate (d('P)dt]. 

A plot of the adaptive loop signal d(¥e)dt from the nominal TAEM case without wind, is 
included in figure H-2 . Notice that over time, d(4 / e)dt approaches zero. 
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Figure H-1 Cross Track Error and Rate Error Algorithm Block Diagram 







155 


NASA 

Naionai At?*0» tfutCS r*nf! 


Report Documentation Page 


1. Report No. 

NASA CR- 187495 


2. Government Accession No. 


3. Recipient's Catalog No. 


4. Title and Subtitle 


Personnel Launch System Autoland 
Development Study 


5. Report Date 
January 1991 


6. Performing Organization Code 


7. Author! s) 

J.A. Bossi 
M.A. Langehough 
J.C. Tollefson 


8. Performing Organization Report No. 


10. Work Unit No. 

906-11-01-01 


9. Performing Organization Name and Address 

Boeing Aerospace & Electronics 
P.O. Box 3999 M/S 82-24 
Seattle, WA 98124-2499 


11. Contract or Grant No. 

NASI- 18762 


12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Langley Research Center 
Hampton, VA 23665-5225 


13. Type of Report and Period Covered 
Contractor Report 


14. Sponsoring Agency Code 


15. Supplementary Notes 

Langley Technical Monitor: 
Final Report - Task 8 


Richard M. Hueschen 


16. Abstract 

The Personnel Launch System ( PLS ) Autoland Development Study 
focused on development of the guidance and control system for 
the approach and landing (A/L) phase and the Terminal Area Energy 
Management phase. In the A/L phase a straight-in trajectory 
profile was developed with a initial high glide slope, a pull-up 
and flare to lower glide slope and the final flare to touchdown. 
The TAEM system consisted of using a heading alignment cone 
spiral profile. The PLS autopilot was developed using integral 
LQG design techniques. The guidance and control design was 
verified using a nonlinear 6DOF simulation. Simulation results 
demonstrated accurate steering during the TAEM phase and adequate 
autoland performance in the presence of wind turbulence and 
wind shear. 


17. Key Words (Suggested by Author(sl) 

Autoland 
Spacecraft 
Control System 
Guidance 
6nnp on 


19. Security Classif. (of this report) 

Unclassified 


18. Distribution Statement 

Unclassified - Unlimited 
Subject Category 08 


20. Security Classif. (of this page) 

21 No. of pages 

22. Price 

Unclassified 

168 

A0 8 


NASA FORM 1625 OCT 86 












